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A Sphere Theorem for the Equations of Elasticity 


By James H. BraMBLE, White Oak, Silver Spring, Maryland, U.S.A.’), 
and College Park, Maryland, U.S.A.?) 


1. Introduction 


In two recent papers [3, 4]%), CoLLins gave expressions for stresses due to 


elastic distributions in an infinite isotropic elastic solid, bounded internally 


by a spherical hollow and a rigid spherical inclusion, respectively. These expres- 


_ sions were given in terms of the corresponding stresses due to the same distri- 
_ butions in an infinite unbounded solid, the singularities of the distributions 
_ lying outside the sphere. 


‘Noo. 


AY oS 


In this paper the more general problem of the sphere with an elastically 
supported boundary is treated. We give here the solutions for the interior 


problem, rather than the exterior problem, although the two are not essentially 


different. The free and fixed boundary cases, given by COLLIns, are limiting 


_ cases of the elastically supported boundary. These are included here also since 
- the forms of the solutions are quite different from those given by CoLtins. The 


~ components of displacement and stress are dealt with in rectangular coordinates 
and because of the symmetry of the sphere the three components of the 
_ displacement vector are given as a single indexed quantity. 


od eae 


2. The Sphere Theorem 


The notation u, and 1;; will be used to denote the components of displace- 


ment and stress, respectively, in rectangular coordinates (%1, Xa, %3). A comma 


preceding a subscript will mean partial differentiation with respect to the 


appropriate independent variable and summation will be implied when an 
index is repeated, e.g. 


E OU; 
i ae ii Os, + Ovqee 


The symbol A will denote the Laplace operator. 


Let u? and 1?; be the displacements and stresses corresponding to a state 


‘of stress in an infinite, unbounded solid with no body forces acting in the 
medium, and such that all singularities lie within a sphere centered at the 


1) U. S. Naval Ordnance Laboratory. 
2) University of Maryland. 
3) Numbers in brackets refer to References, page 6. 
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Now 0; = 6? on v= a, and b° is a regular biharmonic function in 7 > a. Per- 
forming a Kelvin inversion (for biharmonic functions) we obtain the function 
ry b? (a?/r)/a which is a regular biharmonic function in y <a and which is 
equal to 6; on 7 = a. Hence we may say that h, is to be the solution of . 
2 
Alive O00 Sm ree, h,= = 8 (<) on 7 =o. : 
A simple application of our lemma yields (2.4). : 
We now introduce w; as the “regular part”’ of w; in r < a (i.e. ui = uf — u,) : 
and x, t;; as the corresponding stresses, related to w; by Hooke’s law. Then : 
bj = a X; Tj; + &, u;. Using HooKeE’s law we find that : 
| 


Oo 


b= ol gee O' + pu (x; Mi; j + (7%), )| + (a, — @%) U4; , (2.9) 


and that 


ih 1 , , , 
cb, = 4—— O',, where O-= Zz), (2.10) 


Hence 


b= 0, — 00%, = om (x; + (7), ,) + (%y— pay) wy. (2.11) | 


eee 


Now 
' Os | 
X pj) = 2poyr [Ku +r (u)], (2.12) 
or 
r) ? 1 r= 16, 
zy (1% 4) = 2 oy "(5 83) (2.1378 


Hence, (since K = 1/2), integrating (2.13) from 0 to 7, we have 


oe iain © r qa Ve 
Cte nie / o “(p) ao. (2.14) 
6 
Having solved for 7 / in terms of the p,’s we now write (2.11) as . 
4 m) / 1 v | 
(2K—l1)u,+7 57 eee rq Pi = (FU), er I) | 


Integrating (2.15) from 0 to y and using (2.14) we obtain (after a lengthy but 
elementary computation (2.1). | 
Corollary 7. (Fixed Boundary, «, =\0, a, = 1). 
Let u? represent a given state of stress whose singularities lie inside 7 = a, 
Then 


eer tad (Eee) , @ \(1-26)/(3—40) ) 
ioral t= (3—40)7 p (2 h; ;,do, (2.16) | 
i) 
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where 
a 


h,= uf? (=) She 2 — at) [ [os Ale ui(“)] Pye 


“ 0 


represents the state of stress inv < a having precisely the same singularities as 
uy in y < a, and such that uw, = 0 onr = a. We assume that o < 3/2 oro > 3/4. 
This corollary is just a limiting case of the theorem, however, it can be 
derived more simply in a manner analogous to that in the proof of the theorem. 
The details will be omitted. The cases of physical interest are included in the co- 
tollary, however, (2.16) may be extended to all values of o + (3 n—4)/(2n— 3)/2, 
n=0,1,2,..., by analytic continuation in (complex) o. 
Corollary 2. (Free Boundary, a, = 1, %, = 0). 
Let wu? and x; t?; represent a given state of stress all of whose singularities 
lie in < a. If t?; = O(r-*) for large 7, then the formulae 


; ee ee paps eo plat 
; Hats / (2-1) @A).- 22 pifetae, (218) 
7 ’ o ’ 
2 Pica ees Oi,» 
£ (a, — 1) (8-2 ay) (2)" = (a, — 1) B24) (F)" 
a 2B so 
BH ht = 7 a?) [ ll, — ty h;,;iQ do , 
g : (2.19) 
a; mate | 1) |/ot - ; 
z Y 19 (a\2 it) s Q 1/2 2 70 a Mee 0 
i= = 8 (2) are a a) | (2) Al? # (S)| de. Seales BT 
: . (2.20) 


give the state of stress in r < a, having precisely the same singularities as ue 
“inr < a, and such that x,t,; = O0onr= 4a. 

Proof. Under the additional hypothesis, tj = O(r-*) for large 7, Theorem 1 
“is valid for «, = 1, a = 0. For o? < 3/4, a; is complex and Re(a; — 1) > 9. 
_In physically real cases in order to have equilibrium we require that t?; =O(r~*) 
for large rv. Hence, for all cases of physical interest Corollary 2 is valid. The 
_ only values which must be excluded are those for which o = — (n?4++1)/(2n+1), 
ee 2. 3,.... In all other cases the formulae (or continuations thereof) are 
valid. It should be noted that an exterior theorem could have been given as 
well. The formulae in this case can be obtained by changing the limits of the 


“integrals from (0, 7) to (ov, 7). 
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Zusammenfassung 


In dieser Arbeit werden Ausdriicke fiir das Verschiebungsfeld in einer elastischen 
Kugel angegeben, deren Oberflache a) elastisch gelagert, b) eingespannt, c) frei ist. 
Die Verschiebungen werden durch diejenigen Verschiebungen und Spannungen 
ausgedriickt, die in einem unendlichen Kérper unter derselben Verteilung von Sin- 
gularitaten entstehen. 


(Received: Mai 19, 1960.) 


On the Wall Effect Correction of the 
Stokes Drag Formula for Axially Symmetric Bodies Moving 
Inside a Cylindrical Tube’) 


By I-pEE Cuana, Pasadena, California, U. 3, A,%) 


It is well known that the drag on a body moving in an incompressible 
viscous fluid may be determined from Stokes equations for small Reynolds 
numbers. When the flow field is infinite, formulae for computing the drag on 
bodies of various shapes are available [1, 2]*). Recently, a class of solutions of 
Stokes equations for axially symmetric bodies moving in an infinite flow field 
were given by PAYNE and PELL [3]. Due to the nature of the Stokes equations, 
however, the effect of the wall of the containing vessel is generally not negli- 
gible [4]. This is clear from the fact that flow conditions specified at infinity 
often play vital roles to the solutions of Stokes equations. 

The motion of a sphere in a vertical tube filled with stationary fluid was 
first studied by LapENBuRc [5]. Assuming that the radius of the sphere is much 

') Research done under U.S, Air Force Offi 
(638) — 521. 

*) California Institute of Technology, 

8) Numbers in brackets refer to References, page 13, 


ce of Scientific Research Contract Nr. AF 49 
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- smaller than the radius of the tube, he obtained the following formula for the 
_ drag on the moving sphere: 


D=6npaU(1+L->), (1) 


. where a and R are, respectively, the radii of the sphere and the tube; wu is the 


_ viscosity of the fluid; U the velocity of the sphere; and L a constant approxi- 


mately equal to 2-104 [2]. 

The above formula was obtained from an approximate solution of the 
Stokes equations. It gives the drag on the sphere to the order of approximation 
O(a/R). Under the same conditions, we shall show that the drag on any axially 
symmetric body moving inside a tube of radius R is given by the formula: 


2 
D=D,(1+ qa) + (fz) (2) 
where D, is the Stokes drag for the body moving in an infinite domain of fluids, 
and x is a constant approximately equal to 2-203. The second term inside the 
bracket of (2) accounts for the wall effect correction ; for the sphere, Dp=62uUa 
and then (2) reduces to LaDENBURG’s formula (1). 

Our problem consists in finding the drag force on an axially symmetric body 
moving with constant velocity U inside a tube filled with viscous fluid. The 
tube is assumed infinitely long, and the body is moving rectilinearly along the 
centerline of the tube so that an axially symmetric flow field, stationary at 
infinity, is maintained. We choose a cylindrical coordinate, taking the center- 
line of the tube as the z-axis and the radial distance from this axis as the 
radius rv. The origin of the coordinate system is chosen to coincide, at a certain 
moment of time, with a point on the centerline of the body; the exact location 
of this point does not have to be specified. The flow velocity q(7, 2) and the 
pressure #(7, z) are given at this moment by the differential equations [5]: 


uV2q —Vp=0 } (3a) 
divq=0, (3b) 


and the boundary conditions 


= RS ng | at body , (4a) 
q=0 forv=R, (4b) 

and 
q 20 when z—>-+- oo, (4c) 


where i is a unit vector pointing in the direction of the positive z-axis. 
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If the radius of the tube, R, is very large compared with the characteristic 
length a of the body, an approximate solution of (3) which gives the drag 
formula (2) may be obtained in the following manner. We consider the flow 
field in two regions: an imner region which is immediately adjacent to the body 
and an outer region which is far away from the body. The velocity q and the 
pressure # are treated in these two regions by the following successive approxi- 
mations: 

Inner solutions: 


q~ht+h+::-, (Sa) 

b~potpit::. (5b) 
Outer solutions: 

Cie~ By sierra ees (6a) 

b~botbite. (6b) 


The first order inner solution hy and 4, are obtained by letting R > oo. 
The differential equations and the boundary conditions for determining hy and 
py are then 


uV*h, —Vd,=0, (7a) 
div h,= 0, (7b) 
and 
hy=—Ui at body, (8a) 
h, +0 at infinity , (8b) 


and the solution is simply that for a body moving in an infinite flow field, which 
is assumed already obtained. The drag force Dy given by such a solution is also 
assumed known. 

For later reference, the asymptotic expressions for the velocity hy and the 


pressure field 4, valid for large values of r and z, are given below (cf. equation 
(4.9) and (4.11) of reference (Se 


mae (Ea) +O), 0 
By — pte Of) Qh aan et) (9b) 


The above expressions hold true for any shape of bodies with characteristic 
length a. The terms displayed satisfy the Stokes equations and give description 
to the velocity and pressure fields induced by applying a concentrated force of 
strength — Dji on the fluid at the origin. Such solutions are usually referred 
to as the fundamental solutions of the Stokes equations [4] 
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The outer solutions are valid approximations of the exact solutions ina 


_region far from the body; for large tube radius R, this includes the region near 


the wall. In determining the outer solutions using the Stokes equations, the 


_ non-slip condition g = 0 at the wall is used. Near the body the outer solutions 


are not expected to be valid, hence the boundary condition at the body may 
not be used. In the place of the later, we introduce the matching conditions 


_ existing between the inner and the outer solutions. Loosely speaking, this con- 


dition states that the leading terms of the outer solutions when expanded into 
a series for small values of 9 must agree with the leading terms of the inner 
solutions when expanded into a series for large values of 9. When such matching 
is possible, an additional condition is available for the determination of the 
solutions. 

The first order outer solutions g, and f, are thus solutions of the Stokes 
equations which satisfy the boundary condition g, = 0 at the wall and which 
reduce to equation (9) for small values of 9. By direct calculation one can show 
that the proper solutions for gy and #, are given by the equations 


uV2g, — Vb) = Dy d(r) i, (10a) 
div g,=0, (10b) 
and the boundary conditions 
g,=0 fory= Rk, (11a) 
£20 when 2 > 0°, (11b) 


where 6(r) is the 6-function which is zero when r + 0 and is co when r= 0 
and its integral over any volume containing the origin is equal to unity. From 
equation (10a) one sees that the boundary condition at the body (4a) is expres- 
sed by introducing a forcing term Dy 6(r) i in the momentum equation. This is 
equivalent to replacing the body by its retarding force on the fluid. 

The solution of g, may be divided into two parts: gy = gf) + g”). The first 
part, g), consists of terms which are singular at r = 0. These terms arise due 
to the singular forcing function Dy d(r) 7. We shall show by explicit calculation 
that this part is identical to the asymptotic expressions of hy and fy displayed 
in (9). The second part, g®, contains terms which are regular at r= 0 and 
arises due to the non-slip boundary condition gy = 0 at the wall. The second 
flow field is therefore induced by the reflection of the wall. Since the flow is 
axially symmetric. the induced velocity gf) has only the z-component at the 
point r = 0. For small values of 9, one finds (cf. eq. (31) below) 


<) 


ge) = u' i+ O(a). (12) 


\ 
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Next, we consider the second order inner solutions h, and ,. The require- 
ments for h, and #, are then that they should match with gi) which is the 
unmatched part between the first order imner and outer solutions. Thus 


uV7h, —Vp,=0, (13a) 
div, = '0" (13b) 
with the boundary conditions 
hy=0 at body , (14a) 
h,>wi at infinity . (14b) 


The second condition, (14b), is the matching condition mentioned above. By 
considering equations (13) and (14) in conjunction with equation (7) and (8), 
one sees easily that the sum h, + h, is the solution of the Stokes equations 
which correspond to the body moving with a relative velocity U + wu’. The 
effect of the wall, to the order O(a/R), is therefore an apparent increase of the 
velocity of the body in the ratio 1:(1 + w’/U). There is an equal amount of 
increase in the drag force if the momentum integral of the velocity field is 
computed. The drag force on the body in the presence of the wall is then equal to 


DED, (14 5 O( fz). (15) 


This basic formula to be used below was also used by LADENBuRG in obtaining 
(1) based on a different argument. The above formulation follows essentially 
the singular perturbation procedure developed in [6]. In a formal application of 
the theory one should first introduce the small parameter e = a/R and then 
consider the various limits of the flow quantities q and # corresponding to 
different orders of the space variables y and z with respect to e. This systematic 
procedure leads to an inner and an outer expansion similar, respectively, to (5) 
and (6). Such formal development, while leading to justification to many 
arguments given above, is too lengthy to be used in this short paper. 
Now, let us find the value of w’ by considering (10) and (11). We put 


D . 0 
sven [V20 iz Ve (16a) 
Pp=- Dove (16b) 


By inserting these expressions into equations (10), one easily verifies that (10b) 
is identically satisfied and (10a) becomes 


PAVED ee O(F\: (17) 


ce eee 


en elle eee ANP 
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Furthermore, (11a) and (11b) lead to the following conditions for O(r, z): 


(Keo) 0?@ 
a tone == JK. (18a) 
@ +0 when z—>- oo. (18b) 


We define the Fourier transform of a function F(r, z) by the pair 


Fv, a) = / eM Fy, 2) dz (19a) 
and 
oe et Fir, a) d (19b) 


and denote ?® by the function g(r, z). After applying the transform to 
equation (17) there follows 


age 1 dp 


dy? y a 2»? P = d(7) : (20) 
Here 6(7) is a 6-function defined by 
6(r) = d(r) d(x) . (21) 
The solution of equation (20) is 
= 1 A 
yr, A) a aS K,(a r) + Ore I,(A r) . (22) 
The equation for ®(r, A) is then 
2b 1 do ey Ne 
Fri ear BPRO= (ar), (23) 
which gives the result 
2G, 1) =4" Kan) + 3 A ar tant+ gota). (24) 


In (24), A and B are the constants of integration to be determined from 
boundary conditions (18) and Jo, J;, Ko and K, are the modified Bessel func- 
tions. By direct calculation one finds 


TAR) Kg(t R) — 25° [19(A R) Kol R) + 110 R) KAA RY) 
A(AR) = , oR i (25a) 
T(.R) 1,(A R) + > (4A R) — 130 RB) 


WA I-DEE CHANG ZAMP 
and 
AR 
4 
B(A R) = (25b) 


19(R) 1,0 R) + 4* cya R) — 130 RY] 
The z-component of go, u(r, z), is given by the expression 
D 0 
u(r, 2) ==" | lr, 2) — 52 Ot, a)]. (26) 


The Fourier transform of u(r, z) is 


u(r, a) = oe [p(r, A) + 42 Bi, A)] (27) 
= Wr, A) + u(r, 2) , (28) 
where 
= Dy, AY 
wg, A) = — oe [KolA n— Kya ”)| (29) 
and 
u(r, A) = a [A(A%) (Io(A7) + Ar Glan) + BAR) IgA]. (30) 


These two parts of u(r, A) correspond, respectively, to the z-component of 
gi and g?) . The inverse of u(r, 2) is obtained by using (19b): 


/ eon (+ G, o) do; o=AR, (31a) 


== OD 


1 
2aR 


U(r, 2) = 


and w®)(r, z) is regular at the point r= 0; 
: 1 
U(r, z) = 4! 4+ Ol je) ; (31b) 
where w’ is the value of u(r, z) at the point r = 0. One finds 


, 


he woe 
oe rie [p(0, A) + A2@(0, a)] aA 


thee 00 
Sata | MAR) + BAR) gi = 5 Pe [ Ho) ao (32a) 


Sy eae 8 
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where 
co © T,(a) K,(0) — — [I,(c) K,(o) + I,(0) K,(o)] + — 
oe [#0 ee ies ‘i aL 0 2 0 0 1( 1 )] 4 
é ’ Iy(e) Ty(0) + 3 (130) = 1%(0)] (32b) 
=> 2-203 . 
So, finally, 
i uv r * Do 
D Dy (1 + aS eran (33) 


In Table I we list several cases where the value D, is given in reference [3]. 


Table 1 
xed : Drag 
Body Drag (Infinite Domain) (lacide Dube of Radius 2) 

(1) Hemispherical Cup 17-525 wa U 17-525 wa U(1 + 1-956 a/R) 
(2) Flat disc of radius a 16u4aU 16 wa U(1 + 1-786 a/R) 
(3) Sphere of radius a 6a paU 62 ua U(1 + 2-104 a/k) 
(4) Prolate spheroid 8xamauUu 8mamaU (1+ 2-805aa/R) 
(5) Oblate spheroid 8a BuaU ‘8aBhuaU(1-+4 2.805 Ba/R) 


_In the above, a is the radius of the frontal area circle of the body for cases 
(1)-(3). The cross section of the spheroids, (4), (5), intercepted by the 7 — z plane, 
is an ellipse. For the prolate spheroids the foci are at (+ a, 0) and « is equal 
to the value 


1 +1 -1 
ve [-7 (s? + 1) log <*> — s| (34) 


For the oblate spheroids the foci are at (0, + a) and f is equal to the value 
B =[(1—s*) cot-4s + s]-. (35) 


In (34) and (35) s is equal to the ratio of the semi-axis along the z-axis to the 
focal length a, of the ellipse. 
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Zusammenfassung 


In dieser Arbeit wird der Einfluss der Wand auf die Bewegung eines rota- 
tionssymmetrischen Kérpers langs der Achse eines mit zaher Fliissigkeit gefiillten 
Rohres betrachtet. Der Widerstand des K6érpers wird unter der Voraussetzung 
berechnet, dass die Dimensionen des Kérpers gegeniiber dem Radius des Rohres 
klein sind. Wir finden fiir den Widerstand D eines K6rpers, der sich mit der Ge- 
schwindigkeit U langs der Achse einer zylindrischen Réhre mit dem inneren Radius 
R bewegt, Pe 

D=D,[1+ 2. oP cook 
wo D, der Widerstand des K6rpers in einer sonst den ganzen Raum fiillenden 
Fliissigkeit, x eine Konstante etwa gleich 2-203, und yw die Viskositatskonstante ist. 


(Received: June 10,1960.) 


Die Instabilitat der Str6mung zwischen zwei rotierenden 
Zylindern gegentiber Taylor-Wirbeln fiir beliebige Spaltbreiten 


Von Kiaus KirRcHGASSNER!), Freiburg i. Br., Deutschland 


Einleitung 


Diese Arbeit beschaftigt sich mit der zuerst von TAyLor [1]?) untersuchten 
Instabilitat der laminaren, inkompressiblen Str6mung zwischen zwei rotie- 
renden koaxialen Zylindern. Die Stérungen, die von einer gewissen Rotations- 
geschwindigkeit des inneren Zylinders an auftreten kénnen, bilden sich be- 
kanntlich in Form eines in Richtung der Zylindererzeugenden periodischen 
Wirbelmusters aus. 

Um die numerischen Auswertungen mit ertraglichem Aufwand durchfiihren 
zu kénnen, musste TAYLOR, neben anderen Annahmen wie der Kleinheit der 
betrachteten Stérungen, verlangen, dass die Spaltbreite (Differenz der Zylin- 
derradien) klein sei gegeniiber den Kriimmungsradien der Zylinder. 

Die Erweiterung dieser Theorie auf den Fall beliebiger Spaltbreite gelang 
erstmals CHANDRASEKHAR [2] in neuester Zeit. Die von ihm angewandte Me- 
thode zur Lésung des Problems beruht auf einer auch von TAYLOR angesetzten 
Reihenentwicklung der Stérungsamplituden nach Besselfunktionen. 

In der vorliegenden Arbeit wird zurLésung desselben Problems ein anderer 
Weg beschritten, der wesentlich auf der von GORTLER [3] und HAMMERLIN [4] 


1) Aus dem Institut fiir angewandte Mathematik der U 
fiir angewandte Mathematik und Mechanik der DVL an di 
suchung wurde vom Wirtschaftsministerium des Landes B 
*) Die Ziffern in eckigen Klammern verweisen auf d 


niversitat Freiburg und dem Institut 
er Universitat Freiburg. Diese Unter- 
aden/Wiirttemberg geférdert. 

as Literaturverzeichnis, Seite 29. 
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_ benutzten Methode aufgebaut ist. Diese Methode hat gegeniiber der von CHAN- 


_ DRASEKHAR benutzten den Vorteil, dass der numerisch umfangreichste Teil, 


_ ndmilich die Vertafelung der Greenschen Funktionen, von der speziellen Gestalt 
_ der Grundstrémung unabhiangig ist. Die Grundstrémung wird erst im letzten 
Abschnitt der numerischen Auswertungen bendtigt. 


Wie alle auf der Grundlage kleiner Stérungen behandelten Stabilitatspro- 


it " pleme fiihrt auch die vorliegende Untersuchung auf ein Eigenwertproblem, bei 


dem ausschliesslich der kleinste positive Eigenwert physikalische Bedeutung 


_ besitzt. Die Existenz dieses Eigenwertes wird in Strenge nachgewiesen. 


Unsere Ergebnisse konnten in einem Fall mit experimentellen Resultaten 


_ von DonELLy [5] verglichen werden, wobei sich gute Ubereinstimmung ergab. 


1. Das Eigenwertproblem 


Wir untersuchen die inkompressible zahe Str6mung zwischen zwei unend- 
lich langen, um ihre gemeinsame Achse rotierenden Zylindern — ihre Kriim- 
mungsradien bezeichnen wir mit R, und R, (R, < R,) — , wobei wir zweck- 


| miassigerweise das in Figur 1 gezeichnete Koordinatensystem verwenden; u,v, w 


sind die Geschwindigkeitskomponenten in x, y, z-Richtung. 


f 
| 
| 
| 


| 
| 
| 
| 
| 
! 
| 
I 
| 
| 
| 
| 
| 
1 
| 
I 
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Figur 1 
Das benutzte Koordinatensystem. 


Da wir bei den in einer Reihe von Arbeiten wiedergegebenen Ableitungen 
der Grundgleichungen nicht verweilen wollen, sondern nur die benutzten Sym- 
bole und das zu diskutierende Eigenwertproblem anzugeben beabsichtigen, 
gehen wir von vornherein zu dimensionslosen Gréssen iiber. Wir denken uns 
alle von jetzt an auftretenden variablen Langen mit der Spaltbreite h = R, — R,, 
alle variablen Geschwindigkeiten mit der Rotationsgeschwindigkeit des inneren 
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Zylinders FR, w, (w, = Winkelgeschwindigkeit des inneren hip e Os == Wile 
kelgeschwindigkeit des ausseren Zylinders) und den Druck mit o@ Rj oy dimen- 
sionslos gemacht. Dabei steht @ fiir die Dichte. Ferner setzen wir: R;= Ri lh, 
(1) = 1,2). Als weitere dimensionslose Grésse wird in die Gleichungen die 


Reynoldszahl: 
Ry, o,h 
iv, 


Re= 


eingehen. Dabei steht v fiir die kinematische Zahigkeit. 
Die stationare Grundstrémung hat dann, in dimensionsloser Form geschrie- 
ben, die Gestalt (vgl. etwa SCHLICHTING [6]: 


u(y) =A (R,— 9) += (1.1) 
(R, — y) 
mit 
VRE eS (w2/o,) R3 : R, RS (1 2/0) 
Ry thee) 2k, 4 


Thr entspricht eine Druckverteilung #,(y), fiir die sich aus der zweiten Navier- 
Stokesschen Gleichung der folgende Ausdruck ergibt : 


0 rare Oe (1.2) 


Der Grundstrémung werden nun in iiblicher Weise dreidimensionale Sté- 
rungen in Gestalt von in z-Richtung periodischen Wirbelpaaren tiberlagert: 


u(y, 2, t) ) 

(1.3) 
w(y, 2, t) = w,(y) sing z+ e” 
Ply, 2, t) = Polv) + privy) cosa z+ ev! , 


Hierbei gibt o die Zahl der Wirbelpaare in einem Intervall der Lange 22 h und 
y den Anfachungsgrad der Stérungen an. 

Wir beschranken uns im folgenden stets auf neutrale Stérungen (y = 0) und 
bezeichnen fiir jeden Parameterwert o die kleinste Reynoldszahl, fiir die neu- 
trale Stérungen auftreten k6nnen, mit Re,(o). Dann ergibt sich, wie physikalisch 
zu erwarten, dass fiir Re < Reo(o) alle Stérungen mit der Wirbelzahl ¢ gedampft 
(vy < 0) und fiir Re > Re,(o) angefacht (y > 0) sind. Die kleinste aus der Ge- 
samtheit aller Zahlen Reo(o) heisse die hritische Reynoldszahl Ree. 

Geht man mit dem Ansatz (1.3) in die Navier-Stokesschen Gleichungen ein, 
und linearisiert man beziiglich der Stérungsamplituden, so ergibt sich nach 


eae ta 
bad 
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einigen Umformungen das folgende Eigenwertproblem (vgl. etwa WITTING [7]) : 
L{u,(y)] = —2 A a,(y) , 


L4fo,(y)] = 20% 0) (9), (1.4) 


_ Aus physikalischen Griinden interessiert in diesem Zusammenhang lediglich 


der kleinste positive Eigenwert (co) = Rej(c), dessen Existenz wir mit Hilfe 
der im folgenden Paragraphen bewiesenen Hilfssatze nachweisen werden. 
2. Einige Hilfssatze zur Behandlung des Eigenwertproblems 


Bei der Bestimmung des kleinsten positiven Eigenwertes von (1.4) greifen 
wir auf das von GORTLER [3] und HAMMERLIN [4] verwendete Verfahren zurtick, 


' fiihren also (1.4) vermittels der Greenschen Funktionen in 4quivalente Integral- 


_ gleichungen tiber und ermitteln den gesuchten Eigenwert iterativ. Der Existenz- 


nachweis fiir den Eigenwert wird wie in [4] mit Hilfe eines Satzes von 
JENtTzscH [8] gefiihrt, der mit einigen Voraussetzungen iiber den Kern, und 


damit in diesem Fall iiber die Greenschen Funktionen verkniipft ist. Diese 


Eigenschaften sollen im Hilfssatz 1 dieses Paragraphen bewiesen werden. Die 
Aussage des Hilfssatzes 2 benétigen wir spater zur Untersuchung des Eigen- 
wertverhaltens fiir den nur mathematisch interessierenden Grenzfall o > 0, 
das heisst unendlich grosser Wirbelabmessungen. 

Seien G und H die zu untersuchenden Greenschen Funktionen. Sie geniigen 
iiberall im Quadrat Q = {0 < y <1; 0 << 1} mit Ausnahme der Strecke 
y = y den Differentialgleichungen : 


L{G(y, )]=0; L*Aly, 9] =90, (2.1) 


und den homogenen Randbedingungen : 


2 p - ° _ OH(0,9) _ OA(1, 9) 
G(0, 9) = G1, 9) = 9; HO, 9) = HU, 9) =0; SB = 5 = 
(2.2) 


Die Funktion G ist tiberall in Q stetig. Ihre erste Ableitung besitzt fiir y = y 
eine Unstetigkeit. Es gilt: 


Oy oy 


lim { 


e—0 


ZAMP XII/2 


eG 


‘ 
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Die Funktion H ist in Q zweimal stetig differenzierbar. Fiir ihre dritte Ableitung — 
gilt an der Stelle y = y: 


0H) — «, 9) ORY Ae 8,9) |e 2.4 
in See ae ae 


Wir beweisen nun 


Hilfssatz 7: Fiir die Greenschen Funktionen G und H, mit den in den Glei- 
chungen (2.1) bis (2.4) festgelegten Eigenschaften, gelten die folgenden Un- 


gleichungen®) : 


LG eye 0 ine 


Oy oy 
3. H(y,9) <0 ing, 

02H (0, 9) PH(1,9) — 9 = 
4, aor be oy? —< OF Oy? =< ’ y 


Beweis: 1. Die Funktion G kann in Q, ausgenommen auf der Strecke y = y, 
weder ein positives Maximum noch ein negatives Minimum annehmen, da sonst 


die Ungleichungen 
L1G) 2:0 baw 2G) 20 


fiir mindestens einen Punkt in Q gelten. G ist daher in jedem der Teilbereiche 
0<y<y;y Sy <1 monoton. Ware nun G im ersten der beiden Teilbereiche 
negativ, so folgte aus der Monotonie und aus G(0, y) = 0: 


0 fir 0 << ye 


ie 
und wegen der Stetigkeit von G und G(1, #) = 0: 

0G(y, ¥ - 5 

pl 20. fir» By<aae 


Diese beiden Ungleichungen stehen aber im Widerspruch zur Beziehung (2.3), 
woraus sich also G = 0 in Q ergibt. 
2. Langs der Strecken y = 0 und y = 1 gilt daher sicher fiir die einseitigen 
Ableitungen: 
0G(0, 9) 


zen a es 0: 


. aG(1, 9) 
oy — y 


be by 
8) Wie der Beweis zeigt, gilt die Behauptung dieses Satzes auch fiir den Fall, dass der Operator 
L die allgemeine Gestalt 


S0; gel. 


o j d 
L= aly) = + b(y) i + c(y) 


besitzt, wenn nur a, b und c stetig im betrachteten Bereich und a => Ove Oust 
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ee A 


-Angenommen nun, es gebe ein Vy € I so, dass in einer der beiden Ungleichungen 
das Gleichheitszeichen steht, dann folgt aus (2.1) sowie aus der Stetigkeit und 
- Monotonie von G: 


Goi 9,) 20) 0 ys. 1; 


bis he Ran 6 A Si 


A wiederum im Widerspruch zu (2.3). Hieraus ergeben sich aber sofort die Un- 
4 gleichungen der Behauptung 2, und damit G > 0 im Innern von Q. 
3. Setzen wir voriibergehend: 


tee AR ce (2.5) 
so befriedigt K die Differentialgleichung 


L{K] =0; | (2:6) 


a) K ist stetig, und 0K/dy besitzt langs der Strecke y = y eine Unstetigkeit 
der Form (2.3). K kann in jedem der Gebiete 0 < y < yiy<y<lye I nach 
unseren Uberlegungen iiber G weder ein positives Maximum noch ein negatives 
Minimum annehmen und daher in I = {0 < y < 1} héchstens 2 Nullstellen 
besitzen. H selbst verschwindet in keinem Teilintervall von J identisch, was 
sich leicht mit Hilfe der oben erwahnten Eigenschaften von K und Gleichung 
(2.4) beweisen lasst. 

b) Besitzt H fiir ein y = yp) in J m Nullstellen, wobei eine k-fache Nullstelle 
nur 1mal gezahlt wird, so hat K in J mindestens ” + 2 Nullstellen. Um dies 
einzusehen, hat man nur zu beachten, dass H in jedem der durch die Nullstellen 
begrenzten + 1 Teilintervalle von I mindestens cin Extremum annimmt. 
Zwischen je zwei dieser Extrema liegt sicher eine Nullstelle von A, wie man 
sich leicht tiberlegt. Ferner wechselt K zwischen den Randpunkten y = 0 und 
y =1 und den ihnen benachbarten Extremstellen nochmals das Vorzeichen 
(H 0 in der Umgebung dieser Punkte). 

Hieraus folgt aber sofort das Behauptete. 

c) Da aber nach (a) K in J nicht mehr als zwei Nullstellen besitzen kann, 
gilt n = 0. 

Nimmt man nun an, H sei fiir ein y = ¥ im ganzen Interval / positiv, so 
existiert mindestens ein Punkt (yo, Yo), Wo H ein positives Maximum annimmt. 
K ist dann an dieser Stelle negativ, wahrend es auf Grund der Randbedingun- 
gen (2.2) in einer hinreichend kleinen Umgebung der Randpunkte y = 0 und 
y = 1 positiv ist (H 0 in der Umgebung dieser Punkte). Hieraus und mit 
Hilfe der Beziehung (2.3) schliesst man, dass dies nur dann der Fall sein kann, 
wenn K fiir ein y = y, + Yo ein negatives Minimum besitzt, was im Widerspruch 
zu den unter (a) angefiihrten Eigenschaften von K steht. 
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Damit ist gezeigt, dass fiir H in Q die Ungleichung: 


Aly, y) <0 


gilt. 
4. Zum Beweis dieses letzten Teiles der Behauptung nehmen wir an, es 


existiere ein yy € J, fiir welches etwa 


0H(0, 9) _ ¢ 
Oy? = 


ist. Dann gilt aber 0K/dy < 0 in diesem Punkt, wegen H0. Dies ist aber 
wegen der auch fiir K geltenden Beziehung (2.3) und wegen der unter (a) auf- 
gefiihrten Eigenschaften von K unméglich. 

Damit ist Hilfssatz 1 vollstandig bewiesen. 

Wir bezeichnen mit Ly den Operator, der aus L hervorgeht, wenn o = 0 


gesetzt wird: 
ne oe eee & —— ~~. (2.7) 
y fig= y Ol") «(Ry gi}? 


Die Greenschen Funktionen der Operatoren L, und Lj, die den Randbedin- 
gungen (2.2) geniigen, nennen wir G, und H,; dann gilt: 


Hilfssatz 2: Lasst man o gegen Null konvergieren, so gilt gleichmassig im 
Bereich Q: 


Lm |G(y, 9) — Goly, 9) |=0, lim| H(y, 9) = Holy, $) | = 0. 
Beweis: Die Funktion D, = G — Go ist in Q zweimal stetig differenzierbar 
und befriedigt die Differentialgleichung: 


L[D,y(y, y)| = 0? Goly, ¥) . 


Fiithrt man diese Beziehung durch die Greensche Funktion G in Integralform 
liber, so ergibt sich: 
1 
D,(y, 9) = —o* | Gly, 9) Goly, $) dy . 

d 
Bezeichnet Gy das Maximum von Gy im Bereich Q, so folgt, da G nach Hilfs- 
satz 1 in Q positiv ist: 

| Dy(y, ¥) | < 0" Goy ily), (2.8) 
mit 


dy(y) = | Gy, 9) ay. 
H 


or 


VT Ts SF 
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d, geniigt der Differentialgleichung: 
L{ay(y)] = -1, (2.9) 


was man sehr leicht verifiziert. Hieraus gewinnt man aber die Abschatzung: 


a 1 


dy(y) <a, = ofa 1B (2.10) 
2 


Es gilt namlich: 


Da d,(y) fir y = 0 und y = 1 verschwindet, miisste, falls d, > d, an einer Stelle 


gilt, mindestens ein yp existieren, wo d, ein positives Maximum annimmt und 


dy(Yo) > d, ist. Dort bestiinde dann aber die Ungleichung: 
Lid,(y)]< Li[4,]) S—1, 


in Widerspruch zu Gleichung (2.9). 
Aus (2.8) und (2.10) folgt nun: 


| D,(y, y) | <8 


a 
pay 
Da Gp von o nicht abhangt, folgt hieraus: 

lim | Dy(y, #) | =lim| Gy, 9) — Gols, 9) |= 0, 
gleichmassig in Q. 

Den zweiten Teil unserer Behauptung beweisen wir ganz entsprechend. Wir 
setzen: D, = H — H, und erhalten fiir D, die Differentialgleichung : 

LDA 20 Laity) — 0° 414. (2.11) 


Da H, in Q und bei einseitiger Annaherung an den Rand von Q stetig differen- 
zierbar ist, kénnen wir L,[H,] und H, in Q nach oben abschatzen: 


| LlAoly, y)] | =, | Holy, y) < Hoy - 


(2.11) wird wiederum mit Hilfe von H in Integralform iiberfiihrt und | D, | 
nach oben abgeschatzt: 


| Dy, y) | < dy) {2 0? C + o* Hoy} 
mit 


d,(9) = — | Hv. 9) ay 
0 
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d, erfiillt die Differentialgleichung : 


L{d,(y)] =1. 
Setzt man 
L{do(y)] = 9(y) 5 
so gilt 
L{o,(y)] = 1, 


und es ist 0, fiir y = 0 und y = 1 positiv, wie aus der Aussage 4 des Hilfssatzes 1 
folgt. Wendet man die Uberlegungen, die wir oben bei der Abschatzung von d, 
angestellt haben, auf 6, und d, an, so ergibt sich: 


dx(y) < Ea 
a0) (e+ A/RD (a + 1, — D9 


> 
und hieraus fiir Dy: 


: 2 o? 4H 
| Daly, 9) | < sie Bare | 


(oP eT) {0% + Uy (Regiment)*} 


Lasst man o gegen Null gehen, so folgt gleichmassig in Q: 
lim | D,(y, y) | = lim| A(y, ») — Holy, 9) |= 0, 


womit der Satz bewiesen ist. 


3. Berechnung der Greenschen Funktionen und Existenzbeweis 
fiir den gesuchten Eigenwert 


‘Um das in 2 angekiindigte Verfahren zur Bestimmung des kleinsten Eigen- 
werts anwenden zu kénnen, miissen wir die Greenschen Funktionen G und H 
kennen. Zu diesem Zweck ben6tigen wir zunachst ein Fundamentalsystem von 
Lésungen der Operatoren L und L2. Fir L bieten sich hierfiir die modifizierten 


Besselfunktionen 1. Ordnung J , und K, an, deren Wronski-Determinante A, 
den Wert 


4 


A,(r) = — = mit v=o (R, — y) 


besitzt. Hieraus lassen sich aber ohne Schwierigkeit vier linear unabhangige 
Lésungen von L2 bestimmen. Sie lauten: 


il’) = I(r), 


Pal”) = Ky(r) , 


a? emer 


ee a ee 


ty 


f 


= 
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Se 


. 
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mit den Abkiirzungen 


A,(r) = | @ Lilo) Kilo) de, 
Ag(r) = | ets(o)}* do: 


Ag(r) = | ofKx(o)]* ae » 


_ wobei wir fiir 7, den Wert 7, =o (R, — 1) gewahlt haben. Fiir die Wronski- 


Determinante 4 der vier Funktionen 9,(7) gilt: 


. Die Greenschen Funktionen G und H haben dann die folgende Gestalt: 


[Li(72) Ky(7) — Ky (rq) 1,(7)] (11) Kil) — Ky (1) 11(7)] 
fir r<7; 


[14(72) Ky(r) — Ky (72) T,(r)\ (24(1) Ky (7) — Ky (7) 11(7)! 


fiir. 7 =—7;. 


fal?) — 3p- (Aslrs) wal?) — Aulrd) ulA)} 990 


+ {— ul) + gp (Aulrs) al) — Alea) al} alr 
fir x S7; 
Pe) ae 
{72(” a NY [A3("2) ps(”) — Ar("2) pal?)]| y3(7) 
+ [ule + ye LAalr) gale) — Aatea) pale} al?) 


fir 7 = 7, 


ey rer (Rs n= oR — YE =o Rs. 
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Mit Hilfe von G und H lasst sich nun (1.4) in ein System von Integral- 
gleichungen tiberfiihren: 
1 1 - 
" Z PA ey ; xy UY) ‘ 
u(y) = 2 A/ Gly, ¥) "(y) dy, %(y) = 20° al Hy, 9) Baap dy , 
0 0 


oder zusammengefasst 
ni) =f My, 9) (9) ay, (3.3) 
mit dem Kern 


M(y,y) = 402.4 | H(y, 9) 2 Gy, 9) ay. (3.31) 
0 2 


Greift man auf die Aussage des Hilfssatzes 1 des vorigen Paragraphen zu- 
riick, so erkennt man, dass das Vorzeichen des Kerns M durch das Vorzeichen 
von #(y) und von A bestimmt ist; denn H und G besitzen im Quadrat Q keine 
Nullstelle. Wir beschranken uns nun fiir den Existenzbeweis auf solche Grund- 
stromungen, fiir die A < 0 und ~% = 0 gilt. 

Wahrend mit der letzten Ungleichung der fiir die Stabilitatstheorie durchaus 
bedeutungsvolle Fall gegensinnig rotierender Zylinder ausgeschlossen wird und 
damit die Existenzfrage des Eigenwertes hierfiir offen bleibt, kann dagegen der 
Bereich A = 0, das heisst w,/w, = R3/R? ausser acht gelassen werden, ohne 
dass dadurch die Klasse der zugelassenen und fiir die Stabilitatstheorie inter- 
essanten Grundstrémungen eingeschrinkt wird. Man weiss namlich auf Grund 
eines Rayleighschen Stabilitatskriteriums, dass alle moglichen Strémungen der 
Taylorschen Versuchsanordnung, fiir die A > 0 gilt, gegeniiber Stérungen des 
hier betrachteten Typs stabil sind. Dieses Kriterium wurde fiir reibende Fliis- 
sigkeiten erstmals von J. L. SYNGE [9] bewiesen. 

Mit der Aussage von Hilfssatz 1 aus § 2 und den oben gemachten Annahmen: 
uM 20 und A <0 (w/o, < R3/R2) lassen sich jetzt fiir den Kern M aus 
(3.31) sofort die folgenden Beziehungen verifizieren: 


M\y,9)>0 inQ, M(0,¥) = M(1, 9) = M(y, 0) =MVy, 1) =0° 


M ist stetig und beschrankt und daher quadratisch integrabel. Nun verbiirgt 
aber ein Satz von JENTzscH4), dass jeder quadratisch integrable Kern, der 
héchstens in einer Menge vom Mass Null verschwindet und sonst positiv ist, 


4) Der Satz von Jenrzscu wurde zuerst von HAMMERLIN [ 
kleinsten positiven Eigenwertes herangezogen, und zwar fiir de 
Grenzschichtstrémung langs einer konkaven Wand 
schen Zylinder wurde die obige Aussage unter d 
wiesen, 


4] zum Beweis der Existenz eines 
n Nachweis der Instabilitat einer 
gegentiber Gértler-Wirbeln. Im Falle der Taylor- 
er Einschrankung h < R, von Wirtinc [7] be- 
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einen einfachen und vom Betrag kleinsten Eigenwert besitzt. Ferner hat die 
Z zugehérige Eigenfunktion im ganzen Intervall dasselbe Vorzeichen. 

_ Mithin existiert fiir jedes o > 0 ein dem Betrag nach kleinster positiver 
_Eigenwert py = Res von (3.3). v4(y) wechselt im Intervall I das Vorzeichen 
nicht. Dasselbe gilt nach (3.2) und Hilfssatz 1 dann auch fiir u. 


i 


4. Numerische Ergebnisse 


g Zur Berechnung des gesuchten Eigenwertes definieren wir das folgende 
: Iterationsverfahren : 


1 


of"(y) = [ MU, 9) oP) as (4.1) 
0 


Die Folge 
ee] 
poy) 

_ wobei /[v,] ein beliebiges lineares Funktional in 1, bezeichnet, konvergiert dann 
fiir wachsendes » gegen den gesuchten Eigenwert vom kleinsten Betrage, 
vorausgesetzt, dass die nullte Naherung v an der zum Eigenwert gehorigen 
Eigenfunktion beteiligt ist®). Diese Bedingung kann bei der hier vorgenom- 
menen numerischen Durchfithrung des Iterationsverfahrens als erfiillt angese- 
hen werden. 

Da das Verfahren auch dann konvergiert, wenn up das Vorzeichen wechselt, 
die Zylinder also gegensinnig rotieren, sollen sich die numerischen Auswer- 
tungen auch auf diesen Fall erstrecken, obwohl die Existenz des gesuchten 
Eigenwertes hierfiir nicht gesichert wurde. 

Zunachst wurde der Eigenwert 4) = Rep fiir die Parameterwerte 


Rs Ws 


Wy Sma (ony 


= 0) 


- in Abhangigkeit von o berechnet (o = 1, 2, 3, 4, 5). Das Ergebnis ist in Figur 2 
graphisch dargestellt. Der Minimalwert beider Kurven u,, = Re, liegt, wie im 
Taylorschen Grenzfall ( < R;) etwa bei o = 3. Die Wirbeldicke der am frii- 
hesten angefachten Stérungen wachst also proportional zu der Spaltbreite /. 
(o ist die mit / dimensionslos gemachte Wirbelzahl). Wir geben hier noch die 
numerisch ermittelten Werte von Re,,, also die Minimalwerte der in Figur 2 
gezeichneten neutralen Kurven, an: 


Eat 2 (Rep = 1a tihe:*h oe 


\ 


Re,, = 68,30 fir =~ 


5) Zum Begriff der Beteiligung vgl. WIELANDT [10]. 
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Figur 2 
Das Quadrat der Reynoldszahl in Abhangigkeit von der Wirbelzahl oO (Ra = 2h; Ry = 3h; 
(®./@, = 0). 


Ein anschauliches Bild vom Zusammenhang zwischen der kritischen Rey- 
noldszahl Re,, und der relativen Spaltbreite h/R, erhalt man, wenn man, wie 
es in Figur 3 geschehen ist, neben den oben bestimmten Werten fiir Re?, die 
Taylorschen Resultate benutzt, die das asymptotische Verhalten von Re, fiir 
h/R, > 0 wiedergeben. Weiterhin kann die von WITTING in [7] berechnete 
Korrektur des Eigenwertverhaltens von der Ordnung O(h?/R3) verwendet 
werden. Aus Figur 3, in der aus graphischen Griinden Re?, in Abhangigkeit 
von der reziproken relativen Spaltbreite R,/h dargestellt ist, ersieht man nun, 
dass Re;, mit wachsender relativer Spaltbreite fallt; jedoch schwicher als 
linear in R,/h, wie es nach der Taylorschen Theorie zu erwarten gewesen ware. 

Die kritische Re-Zahl wurde ausserdem fiir verschiedene Verhaltnisse der 
Winkelgeschwindigkeiten @®2/@, bestimmt. Das Minimum von Re; liegt, wie 
Figur 4 zeigt, in beiden Fallen in der Nahe des Wertes w,/w, = — 0,1. Die 
Grundstrémung ist daher, fiir die betrachteten endlichen Spaltbreiten im 
Gegensatz zum Taylorschen Grenzfall h > Ry, bei variablem (®»/@, am friithe- 


sten instabil, wenn der dussere Zylinder mit der Winkelgeschwindigkeit 
®, = — 0,1 a, rotiert. 


- chungen an Stelle der hier ermittelten Reyno 
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Figur 3 


_ Das Quadrat der kritischen Reynoldszahl in Abhangigkeit von der reziproken relativen Spaltbreite 


(o = 3; w,/@, = 9). 


Die numerischen Ergebnisse erlauben fiir den Fall Ry = 2 h einen Vergleich 
en Resultat von DonELLY [5], der in seinen Untersu- 


mit einem experimentell 
ldszahl die kritische Taylorzahl Ta 


_angibt. Diese ist definiert durch die Gleichung 


[1 — Ri/Re? 


ated (2 ee is (1 = (9/4) (1 — R3 Og) RY 4) 


v 
Fiir @,/@, = 0 und h = R,/2 ergibt sich hieraus: 


ae 64 / Rim, \2 
Fas (") . 
Umgerechnet auf die Reynoldsche Zahl Re = k, o, hiv lautet DoNELLY’s Er- 


gebnis in diesem Fall: 


Chie)? .... m8 
Tay, = 5 ( ; ye S Ret, = 32780 + 130. 


Durch Multiplikation der oben berechneten Grosse Re;, = (68,30)? mit dem 
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Figur 4 
Das Quadrat der kritischen Reynoldszahl in Abhangigkeit vom Verhaltnis der Winkelgeschwindig- 
keiten w,/w@, (6 = 3; Ry = 2h; Ry = 3h). 


Faktor 64/9 ergibt sich: 
64 2 
<r Re, = Ta,, = 33173. 


Das errechnete und das experimentell bestimmte Ta,, weichen also um etwa 
1° voneinander ab. 

CHANDRASEKHAR gibt in [2] als berechneten Wert Ta kr = 33100 an, dessen 
prozentuale Abweichung vom experimentellen Wert um ein Geringes kleiner 
als der hier ermittelte Wert ist. Letzterer liesse sich jedoch durch Verwendung 
einer genaueren Integrationsmethode beim Iterationsverfahren weiter verbes- 
sern. 

Zum Schluss soll hier noch das Verhalten des Eigenwertes vy) = Re? fiir 
a > 0 bestimmt werden, das heisst fiir den nur mathematisch interessierenden 
Grenzfall unendlich gross werdender Wirbelabmessungen. 

Wir betrachten zunachst die Greenschen Funktionen in diesem Grenzfall. 
Dazu ziehen wir dir Aussage des Hilfssatzes 2 heran, wonach die Funktionen G 
und H fiir o + 0 gleichmissig gegen die Greenschen Funktionen G, und A, 
der Operatoren Ly und L? streben, Fasst man nun o* uy in (3.3) zu einem neuen 


a 
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Eigenwert Mo Zusammen, so besitzt w, im Grenzfall o = 0 einen endlichen 
-Grenzwert. Dies folgt aus der Tatsache, dass sich der Eigenwert bei gleichmassig 
‘Stetiger Veradnderung des Kerns stetig andert (vgl. BUCKNER [11]). Der Wert 
“von fy fiir o = 0 ergibt sich mit Hilfe des oben beschriebenen Iterationsver- 
ahrens aus der Integralgleichung (3.3), wenn man dort G durch Gy und H 


durch H,§) ersetzt. 


* 


oe 


ae 


Bp Fiir das asymptotische Verhalten von uy selbst folgt dann: 


Mo=— fir o>0, (u=const), 


‘ein Resultat, das insbesondere auch fiir den Taylorschen Grenzfall h < R, 


4 


gilt und iibrigens auch fiir die Gértlersche Theorie der Grenzschichtinstabilitat 
_an einer konkaven Wand. Im ersten Fall wurde es von WITTING [7], im zweiten 


_Falle von HAMMERLIN [12] bewiesen. 
; 
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6) G, und H, besitzen eine ahnliche Gestalt wie G und H. Sie gehen aus diesen hervor, wenn 
-man in (3.1) 7, durch Ip, K, durch Kg und 7 durch 9 = R, — y ersetzt. 
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Summary 


An investigation is made of the instability of stationary flow between two 
rotating cylinders for any gap distance. Previous work on this problem, viz. the 
Taylor theory, has been done with the restriction of small gaps. Moreover, in 
contrast to the recent work of S. CHANDRASEKHAR [2], who has also analysed this 
problem, the method used is independent of the basic flow and is therefore valid for 
all flows between two coaxial cylinders. In fact, it can even be used for the analysis 
of the flow in a curved channel. 

The major item of physical interest is the critical Reynoldsnumber, i. e. the 
value of the Reynoldsnumber where small disturbances are amplified for the first 
time. The value of this parameter is determined by the smallest positive eigenvalue 
of the boundary value problem. An existence proof is made for this eigenvalue for 
any wavenumber and for all possible cases of cylindrical flow, with the exception 
of the case where the cylinders rotate in opposite directions. The results are depicted 
in a convenient form, where the critical Reynoldsnumber is the dependent variable 
and the gap distance and the angular velocity ratio of the two cylinders are the 
independent variables respectively. 

A comparison with experiment is made for the case where the ratio of the two 
radii is 2:1 and the outer cylinder is at rest. The agreement with the theory is 
good; the noticeable error being approximately 1%. 


(Eingegangen: 20. April 1960.) 


Photoelastic Calculations by a Complex Variable Method 
By Joun A. Lewis and Henry O. Potrak, Murray Hill. N. J., USA) 


1. Introduction 


From photoelastic observations of isochromatics and isoclinics, one may 
determine the principal stress difference, p — q, and the principal angle « at 
each point of a slab of transparent, elastic material in a state of plane stress. 
The shear stress o,, and the difference of the normal stresses o,,, 6,, May be 
calculated immediately from the relations ve 


20,,= (b—4q)sin2a, (1.1) 


‘ On. Oyy re (p aa q) cos 2 (1.2) 
or, In complex form, 
£0 =0,,—0,, + 210,, = — (be g) ee, (1.3) 


where w is the ‘conjugate stress deviator’ (see, e. g. [1]) 2). 


1) Bell Telephone Laboratories, Inc. 
*) Numbers in brackets refer to References, page 37. 


. 
| 


“ 
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z In order to determine the normal stresses themselves, it is necessary to 


_ make further experimental observations (e. g. of the thickness change of the 
_ loaded slab) or to perform further calculations, using, for example, the given 
_ values of # — q and « and the values of the stresses at one point. It is the 


- latter problem with which we are presently concerned. 


TOL SS 


e 2. Shear-Difference and Relaxation Methods 


The equilibrium equations for an elastic body in a state of plane stress are 


ay (2.1) 


: ame oy 
: WOny OCyy _ 
q Ox + Oy waa , (2. 2) 


or, in terms of the conjugate stress deviator w and the mean normal stress 
oO =(0,, ee 0,,)/2, 


Oo ow 
Oz’ (2.3) 
where 
Leta ke game ke Chas 67/2 ere) 
Oz de Sv, ae are <a Ss +4 55): 


For compatibility o must be a harmonic function. Thus o has the form 
o = 2 Re @(z) = D(z) + D(z), (2.4) 
where @ is an analytic function. Equation (2.3) then becomes 


B= 2. (2.5) 


If w is given, Equation (2.3) or (2.5) shows that o or ® may be calculated up 
to an integration constant. Equation (2.5) yields 


Ble) = Dla) + [Ge ae (2.6) 


“0 


This is the technique employed (in real form) in the so-called shear-difference 
method. It requires the numerical differentiation of w, which is given experi- 
mentally. This may place undue emphasis on random experimental errors when 
the mesh spacing in a finite difference scheme is small. 

This source of error may be avoided by calculating the harmonic function 
o in the region directly from its boundary values, e. g. by relaxation. The 
boundary values may be determined from the surface tractions (X, Y) and the 
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stress deviator w from the formula 


idF =i(X —i1Y)ds=odz+wdz, (2.7) 


for the x- and y-components of surface traction on the line elements ds of the 
arc z = 2(s). This method of ‘harmonization’, however, requires detailed know- 
ledge of the surface tractions, information which may not be available. 

The shear-difference method could be used to determine boundary values 
of o for relaxation in a subregion, but one would still be limited in accuracy by 


the numerical differentiation required. 
In the following sections we shall obtain formulae for ® which avoid this 


difficulty by the simple application of the Cauchy integral formula on circles 
in the region. 


3. The Single-Circle Method 
We first note that Equation (2.5) implies that w has the form 
w= 2@"(z) + Y(2z) , (3.1) 
where Y is an analytic function. This may also be written as 
(2 — 29) O'(z) + Pz) + 2, D'(z) = w/z) . (3.2) 


Now consider the values of the above functions on the circle C: |  — z| = R, 
where C and its interior lie within the region D under consideration. Since on C 


oe ie 
eS 4 0iaa oe ) 
Equation (3.2) may be written on C as 
R? (6 — zo)? B'(C) + WC) + 2% B'(C) = w(C) - (3.3) 


We observe that both Y(z) and z)®’(z) are analytic functions on C and its 
interior and hence contribute nothing to a contour integral of Equation (3.3) 
around C. Furthermore R is a constant, so that we find on integration, using 
the Cauchy integral formula, 


Rf @(2) ac ; a ore 
zn | Fa HP O'e) =z : w(C) at , (3.4) 
Or 
1 27 
®' (2) = = | w(% + Re) 9 do, (3.5) 
é 
where 


CS2+ Re? 


4 
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_ Parenthetically, note that, if we multiply Equation (3.3) by (¢ — 2)" and 
integrate, we obtain the identities 


Wea 


22 
3 Je (athe an! dO =O. Tor 9 = 2.3. ies (3.6) 


0 


_ These might serve as checks on the consistency of the experimental data. It is 
also interesting to note that Equation (3.5) may be obtained by integrating 
_ Equation (2.7) around C. Equation (3.5) thus merely states that the resultant 
_ force on C must vanish. 
j @’(z) can thus be calculated throughout the interior of D by laying down 
circles of convenient size around every point at which its value is desired. The 
_ function @(z) may then be determined, up to an integration constant, by 
taking a line integral along any path in the region, so that 


@(z) =, + i @' (2) dz, (3.7) 


where the integration constant ®, might be the value of @ at the origin, taken 
to lie in the region. The real part of ®,, which is all which is required, may be 
determined by inspection if the stresses are given at some interior point. 
Usually, however, the point at which the stresses are known lies on the boundary 
of the region (e. g. on a boundary free of surface tractions). In this case Re Dy 
might be determined by extrapolation of ® — @, to the boundary and com- 
parison with the given value at the point. 

This method uses the given data rather inefficiently, a set of values of w 
on a given circle being used for the determination of the value of @’ at but a 
single point. We shall refer to this method as the ‘single-circle’ method. In 
the next section we shall develop methods which use the data more efficiently 
at the expense of slight additional complexity. 


4. The Two-Circle Method 


Equation (3.3) can also be written in the form 
[R? + 2 (C — 2)] (CE) + E — 20) P(E) = (F — 20) WUC) - (4.1) 


If we multiply this equation by (¢ — 2)", where z lies within C, and integrate 


around C, there results 
; 1 — 2) w(¢) ae 
Re + % (¢ — 2) 92) + @—m) He) = zr; f SRP 4 


C 


by the Cauchy integral formula. This relation holds on the interior of any circle 


ZAMP XII/3 
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C,:|¢— 4%, ) = R, in the region D. In particular, if we take two circles Cj, C2, 
having some portion of their interiors in common, we obtain two independent 
equations for @’(z) and ¥(z) for z in this common region, viz, 


[Re +% @—a)lO'@) + (e—4)¥@)=W, k=1,2, (4.3) 
where : (c ) w(t) at 
— 2) w 
“= 271 /f (=—2 ; eae 
Ck 


These equations, of course, are not solvable for values of z at which they are 
not independent. For example, if the circles are concentric (z, = z, = 2), at 
the center z = 2) Equation (4.4) reduces to Equation (3.4). In this special case 
the equations obviously cannot be solved for Y/(z9), since its coefficient vanishes. 
In general, the equations are independent if the coefficient determinant 


Ri 2, (2 — 21), 2= % 


A(z) = (4.5) 


RE + 2 ( — 22), z— 2 | 


is non-vanishing. It turns out that if one circle lies entirely within the other, 
A(z) always has one zero within the smaller circle, i. e. in the region of appli- 
cability of Equation (4.3), and one outside, both of them lying on the line 
passing through z, and z,. If the two circles touch, A(z) has a double zero at 
the point of tangency and finally, if they intersect, the zeros lie at the inter- 
section points. 


The use of concentric circles C,: | — z | = Ry Ca: [C= & |= Ry, wit 
R, < R,, leads to an especially simple expression for @'(z) namely, 
, W, — W, 
cD) (z) ==y RR ) (4.6) 


for | z— z| < R,. Obviously this formula also holds at z= 2) by Equation 
(3.5). It yields values of @’(z) all over the interior of C ;, We may now compute 
@(z), up to an integration constant, by evaluating a line integral on @'(z), given 


by Equation (4.6) over any path lying within C,. If we integrate Equation 
(4.6) there results 


Ou 
D(z) = D(,) + 2 4.7 
R3 — R? 
with 
1 es ea = = 
Q -/ Watites e / (C — x) w(t) logs “#8 de . (4.8) 


We have now obtained a rather simple formula for © itself which holds all 
over the interior of the circle C,: | z— Zy| < R,. Note that for all z in C,, 


7 


ep 
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- one uses the same sets of values of w on C, and C,. Once w has been evaluated 


_ by experimental observation on a set of suitably chosen circle-pairs, covering 
_ the region, ® may be computed directly from the above. If we set 


C= 2, + RK, e'” ween ©, and ©, Bae Ye? 0 reel, 


4 gh Ry’ w,(0) = w (zy + R,e°), 
we obtain 
27 
2 
Q,(r, 9) = Fh [2 w(6) log (1 — 7 e-") a0, (4.9) 
6 
D,(7, ¢).= ae je w2(6) log(1—77 et (P-9)) dQ, (4.10) 


Similar, but considerably more intricate, formulae are obtained when the 
_ circles are non-concentric. Although for a region of some particular shape it 
_ might be advantageous to use such a configuration, in general these formulae 
_ do not appear to be particularly useful, especially when one recalls that some 


neighbourhood of the zero of A(z) must be excluded. 


5. Determination of the Integration Constant 


We now assume that the computation implied by Equation (4.7) has been 
carried out over a suitable set of covering circle-pairs in the region. If the 
circles overlap, the set of integration constants, one for each circle-pair, can 
obviously be reduced to a single one, ®) say. We assume that this ah) been 
done. We now have a function 


P(2) = D(z) — D,, 


given all over the region D and we wish to determine ®, or, more precisely, its 


real part. 
As has been mentioned, it is necessary to add other information than the 


_ photoelastic data, given by w, to determine the stresses completely, i. e. in the 


present case to determine Re ®,. Usually, in practice one knows the total load 
on some segment of the boundary. For example, on any segment of a free 
boundary the total load is zero. We now develop a method for the determination 
of Re@®, which depends on the knowledge of the resultant force on some 
finite segment of the boundary. 

By Equation (2.9), the resultant force on an arc AB is 


B B 
Vee 2 | Re B(2) dz al w(z) dz 
A A 
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or, substituting ® = @ + @, and solving for Dy, 
B B 
2(B— A) Re@y =i F,, — 2 | Re D2) dz + | w(2 dz. (5.1) 
A A 


If we take A and B to be two points on the boundary between which the 


resultant force Fy, is known, Equation (5.1) yield Re ®p, since both ®(z) and 
w(z) are known in the region. This completes the solution of the problem. 


6. Further Observations on the Method 


The derivations given in Sections 3 and 4 depend for their success on the 
fact that on the periphery of a circle, z is a very simple function of z. If 2 is 
the center and R the radius of the circle in question, then, on the periphery, 

2 
peje (6.1) 


Z— 2 


Using this relation, we are able to do contour integrals of expressions involving 
z by the calculus of residues, and this is the key to our technique. Now on any 
analytic curve, z may be considered as a function of z, which function will be 
analytic in a neighbourhood of the curve (although it represents 2, of course, 
only on the curve). If on any curve y, z should happen to be a rational function 
of z, we may again obtain ®(z) from Equation (3.1) by contour integration. 

Unfortunately, it turns out that we have already exhausted what can be 
done in just this fashion. For we have the following: 

Theorem. Let z = S(z) on a closed analytic curve y. Then S(z) is a rational 
function of z if and only if y is a circle, in which case S(z) is given by Equation 
(6.1). 

This theorem appears in [2], which paper also studies other properties of the 
function S(z) for arbitrary curves y. The most important for our present 
thinking are that S(z) must always have singularities inside the closed contour y 
and that S(z) cannot have a single-valued integral. 


The next simplest case after circles which one can consider is that of 
ellipses. On the ellipse E: 


(% — %)? 


ae +e se DL, (6.2) 


z is given by the following function of z: 


S(z) = % + (z — 2) Lee SeeM = @ = Pie Ho)" 
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Hence we may write, after some manipulation, 


i @’(C) Zoe PC) dt + b? + a? i @’(¢) dt 
C= 2 | 


(€—2) (C—4) “° * Ba 


E 


2ab_ * V1 — (a? — b%)/(E — 2)? dt ( w(t) dé 

ise) aii (Cane 7) 
If we now take three confocal ellipses, then each of the three integrals on the 
left-hand side of Equation (6.3) will be the same no matter which ellipse we 
are on. This is true since z, the center of the ellipses, and a? — 6?, up to a 
possible minus sign the square of the distance between the foci, will be the 


same. Hence three equations of the form Equation (6.3) can be solved simul- 


taneously for, in particular, the second integral, which equals 21 @'(z). We 
can do this if the determinant 


1 B+a 2a, 


1 b2+a2 2ayb, 


1 B2+a2 245d; 


where a; = bF — c2, does not vanish, and this is true in general. 
Thus we may use three confocal ellipses Z; and obtain @(z) up to a constant 
@(z,) at any point interior to all three by forming a linear combination of 


i w(C) riz 
2m ‘i ae S log, & 
Ej 


as was done in Section 4 for two circles. The advantage which the three 
ellipses might offer in some cases over two circles is that a larger portion of a 
long thin region:may well be within three ellipses than within even several 
pairs of concentric circles. 
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Zusammenfassung 

Die komplexe Darstellung der allgemeinen Lésung fiir den ebenen elastischen 
Spannungszustand fiihrt zu einer einfachen Formel fiir den Mittelwert der Normal- 
spannung im Priifstiick, die sich spannungsoptischer Messungen auf zwei Kreisen 
bedient. Diese Formel vermeidet die Messfehler betonende numerische Differen- 
tiation und bendtigt die Kenntnis der Randspannungen an jedem Randpunkte 
nicht. 
(Received: Mai 5, 1960.) 
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Design of Large Permanent Magnets with Rotationally 
Symmetrical Poles’) 


By DimirrijE UGRIn-SpaRAc, ‘Zagreb?), Yugoslavia 


1. Introduction 


All difficulties in design calculations of a large permanent magnet arise not 
only from nonlinearity of the demagnetization curves, but also from the fact 
that permeance of the air is not negligible as compared to that of the magnetic 
material. 

To be precise, we must distinguish between two kinds of losses: fringing in 
the gap and leakage from the surfaces of the magnetic material and polecaps. 

Classical methods of design calculation use HopKINSON’s leakage factor, 
obtained by conformal mapping for appropriate configurations. Though the 
leakage factor method is a very simple one, and accurate within practical limits 
(see for instance references [1]) and [2]) it does not give any deeper insight 
into the refined geometrical problems of large permanent magnets. 

Another method, that can be called the distributed parameter method, will 
be developed and applied here. It concerns above all the design of large perma- 
nent magnets for studying nuclear magnetic resonance, but it is also applicable 
equally well to permanent magnets of any size, provided that initial conditions 
(as defined in Section 5) are in the proper range. 


2. One Approach to the Problem 


Consider the pair of identical rotationally symmetrical magnetic poles with 
axial cross-section as indicated in Figure 1. Making differential shifts from 
radial cross-sections 1-2 and 3-4 to radial cross-sections 5-6 and 7-8 respec- 
tively, increase in magnetic flux of either pole will be 


d® =U aP (1) 


where U is the Magnetic potential difference between rings 1-2-5-6 and 
3-4-7-8, dP permeance of the corresponding air shell. For the sake of simpli- 


') This work was resulted from a preliminary design calculation of the large permanent magnet 
for high resolution nuclear magnetic resonance, that will be constructed for the Department of 
Structural and Inorganic Chemistry, Institute « Rudjer Bo&skovié », Zagreb. 

*) Institute « Rudjer Bo&kovié », Zagreb. 

’) Number in brackets refer to References, page 52. 
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_ city in making an approximative design calculation we shall boldly assume 
_ that the lines of force 1-3, 24, etc. are parts of circles with centers in the 
_ origin O. This idea was borrowed from [3], p. 107. It should be noted at once 
that to assume circular lines of force leads, at least in the vicinity of the gap, 
to an overestimation of the size of the magnet, but the usefulness of such an 
approximation will be further justified. 


/ ¢ 7 e =-4 ~ 


Ta ae CR No laa _| slope ky 


slope ky *s. 
SA ooh bs 
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Figure 1 
Axial cross-section of the rotationally symmetrical magnetic poles, with indicated substitution 


(dot-dashed lines) for curved tapering of the polepieces. The figure also explains symbols used in 
the text. 


Modifying the expression for the permeance of a spherical air shell, derived 
from [3] (page 109, Equation (13)), we can obtain 


Tt fo d(x? + 7°) 
arsinh (4/7) (2) 


aP = 


The second equation that we need for the design calculation of the permanent 
magnet would read 
aU =—2H,, dx, (3) 


where H,, is the magnetic field intensity in the magnet. 
Differential equations (1) and (3) could be rewritten as follows 


a® aP 
ea ae a 
oO) ew (3a) 


ties 
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These equations together with the demagnetization curve of the selected 
magnetic material entirely describe the situation that we encounter in designing 
large permanent magnets with rotationally symmetrical poles. 


3. On the Demagnetization Curve for Large Blocks of the 
Anisotropic Magnetic Material 


ANDREW and RusHwoRTH [1] have already shown that there are great 
discrepancies between the demagnetization curves of small and large aniso- 
tropic magnetic material specimens. They found decrease in the average BH- 
product for their Alcomax III poles to be 1.5 m.g.o., namely from 4.75 m.g.o. 
for small specimens this product has fallen to 3.25 m.g.o. for large blocks. 
Hence, when dealing with anisotropic magnetic material, we shall generally 
reduce the factory data for (BH),,,,, by a factor of 0.685. 

Another large permanent magnet (University of Leeds magnet) was made 
from Columax magnetic material. A factory leaflet gives the range 600 to 
640 Oe for magnetic field intensity at the (BH),,,, point. Now, from induction 
in the gap and dimensions of this magnet we find 0.8 to 0.85 to be the range of 
the factor, by which we must reduce the magnetic field intensity at (BH)... 
point (and anywhere along the demagnetization curve). It is perhaps a matter 
of taste which number will be chosen from this range. We shall take 0.85. Thus, 
decrease in the magnetic induction would be determined by the factor 
0.685/0.85 = 0.8. 
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Reduced demagnetization curves for large blocks of the denoted magnetic materials 
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Figure 2 shows Alcomax III and Columax reduced demagnetization curves, 
prepared for future reference. 

i From the above discussion we can conclude that magnetic material data, 
given for small specimens of anisotropic material, may be very imprecise for 
‘large blocks, making any method of design calculation only approximative. 
This fact presents the first justification of the approximation introduced in 
the Section 2. ~ 


- 
a 
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3 
: 4. Distributed Parameter Method of Design Calculation for Large 
4 Permanent Magnets With Rotationally Symmetrical Poles 


- 4.1. Calculation of the polecaps. One use of the polecaps is to match induc- 
tion B,, in the interior of the polepieces with induction B, within the gap. 
_ Another use of the polecaps is to homogenize the field within the gap, as is the 
case, for instance, with large magnets for nuclear magnetic resonance. 

We shall investigate only this matching property of the polecaps. The 
"quotient 
m= —* (4) 


_ will be named the induction ratio. Supposing polecaps with infinite permeance, 
the magnetic potential difference between the surfaces of these polecaps will be 


U,=2dyH, = 2dy—*, (5) 
Fo 

where 2 d, is the gap length according to Figure 1, and H, the magnetic field 
_ intensity within the gap. To obtain normalized formulae, let the half gap length 
dy be the length unit. Referring to Figure 1, we shall use the following designa- 


tions: 


. a d l é Y 

; &= —_; (i et hea ————. 735 | 

_ nae hy d d d 
) ) 0 0 0 | (6) 
Kier lr eee Cr = Cor: op em Sy | 


ki = 


I 


ee a yg? REG, Re ae 


Total flux, flowing from either polepiece into the attached polecap, divides into 
two parts: flux ®, that flows through the gap, and flux ®, representing the 


leakage flux from the polecap lateral surface. 
As it is made in reference [3], p. 68, fringing in the gap can be taken into 
account enlarging diameter of the gap by one gap length 


®, = B, (% + dy)? x . 
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This is valid for cylindrical polecaps if the gap ratio 9, = 5, but for conical 
ones with slope k, about 1 we shall extend this range to 0) = 3. 
Further, leakage flux ®, can be found by means of equations (1), (2) and (5): 


dy ae 
0, = u0 0, VETO) 2 apn BT, 
é, 


arsinh (4/7) 


where 


6 
ky + &/o 
[= fe dé . 
a /1 + (€/e)? arsinh (E/g) ; 


Denoting 7 = é/o, one part of the integrand would be 


i] 
% = ——— : 
V1 + 7? arsinh y ”) 


The graph of this function, which will be important later, is shown in Figure 3. 
It can be proved that 7 < 6/0, if ky < gy. This latter condition will be fulfilled 
for every practical case. Thus we may expect 7» to be less than, say, 7/11 = 0.64. 
A first approximation for the function x = x(n) for the given range of 7 can be 
seen from Figure 3: x = 1. Approximation for the remaining part of the inte- 
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Figure 3 
Plot of the function % = (7). 


grand is then Roo/é, where 9 = gy + hy (5 — 1). If this bold simplification is 
allowed, integration is then entirely simple: 


I = (AB + 1) (6 1) = Ry (0p — Bho. (8) 


We should note that, as the inte 


grand is enlarged, values calculated from (8) 
will be increased. 


, 
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: 

q Induction on the boundary between the polepiece and the attached polecap 
“was our final aim: 

A ®,+ @ (eget) OT: 

-4 B = 1 27a 0 

3 me Paley + hy (O— 1) 


a 
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and with the definition equation (4) 


(Qo + Ao (6 — 1))? 


CORES eee ” 
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_ Since the approximative value for J was enlarged, the approximative value 
for m is diminished. Hence the use of the expression (9) for m leads to additional 
oversizing of the magnet. 

Indeed, although it is not here proved, it may be hoped, that derived 
formulae cause oversizing of the magnet at least as long as ky S Qo. 

_ Taking the gap ratio oy as constant and slope k, as parameter, the induction 
ratio m from (9) can be plotted against 6. Two families of curves, corresponding 


to 0 = 3.7 (continuous curves) and g) = 7 (broken-line curves) are shown in 
Figure 4. 


Figure 4 


Families of curves corresponding to the function m = m(Qp, kp, 0) with ky as parameter and gap 
ratio @) = 3.7 (full curves, left and lower scales) and @) = 7 (dashed curves, right and upper scales). 


If one family of curves is intersected by a straight line m = constant, and 
0, calculated for pairs of %y and 6 thus obtained by means of the equation 
01 = Qo + Ro (6 — 1), one can construct another family 0, = 0,(m, 6), where m 
now plays the role of a parameter. Such new families are drawn in Figure 5 
for values 0) = 3.7 (continuous) and @) = 7 (broken line). What is characteristic 


I 
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and interesting for curves of these families is the existence of the ae 
value of 0, at the fixed induction ratio m. Since, however, this fact is base 
upon the rather dubious assumption of circular lines of force, we shall try to 
give a physical explanation. 


Figure 5 


Families of curves corresponding to the function 01 = 0:(Qo, , 0) with m as parameter and gap 
ratio @) = 3.7 (full curves, left and lower scales) and @) = 7 (dashed curves, right and upper scales) 


Imagine Figure 1 representing one limiting case: 7, > 75, d, = dy. Let d, be 
slightly increased with 7, and 7, unchanged. Increase in the cone lateral surface 
will be negligible by comparison with the increase in the length of the lines of 
force. Hence, leakage flux will be decreased and the induction ratio m enlarged. 
To retain the same induction ratio, we must diminish 7v,. On the other hand, let 
%; > 7%, but d, much greater than dy. Decreasing d, slightly, with 7) and 7, 
constant, lines of force lengths will be decreased negligibly as compared with 
the decrease in the cone lateral surface, thus diminishing the leakage flux and 
enlarging the induction ratio m, Again, to maintain the same induction ratio, 
we must diminish 7,. 

Minimal values for 9,, with m and Q fixed, do not necessarily cause minimal 
volume of the polepiece, as will be shown. But deviations from minimal 0, being 
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unimportant in practice, we shall use the functions 6 =6(m, 09) and ¢,=0(™, 09), 

defined by minimal values of 0,, taking 9) as parameter. In this manner, from 

_the families of curves drawn in Figure 5 we arrive at the families of curves 
shown in Figure 6. 

: 10 


4] 


Figure 6 
Families of curves corresponding to the function 6 = 0(m, Qo) (full curves, left and lower scales) and 
0; = 01(m, Qo) (dashed curves, right and lower scales), defined by minimal values of @,, with gap 
ratio 0g as parameter. 


4.2. Calculation of the polepieces. Leaving a detailed discussion of the optimal 
design for the next Section, we shall deal here with polepieces having a fixed 
operating point, i. e. where throughout the magnetic material induction B,, is 

“constant. To compensate for leakage flux, radial cross-section of the polepiece 
must increase from the gap-end towards the yoke-end. As a first attempt to 
find the shape and dimensions of a constant induction polepiece, we shall 
neglect the influence of the yoke upon the field distribution. Axial cross-section 

’ through the field and the poles is therefore as shown in Figure 1. Let us sup- 
pose: 1) polecaps with infinite permeance, 2) an iron yoke with relative 

_ permeability m,, (from 500 to 1000), and 3) the validity of the equations (1a) 
and (3a) for this case. 

To determine the length of the polepiece, we must use the well known 

- magnetic circuit equation 


Qipw, 21H, 2b, =, (10) 
- where the new symbols have the following meanings: 


L_ half length of a ‘C’ arm of the double PG. yoke, 
H, magnetic field intensity inside the yoke. 
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Analogously to the equations (4) and (6) we shall define similar nondimen- 
sional quantities 


N= 


Hy B, (G) : Te Bin E 11). 
He (ey an a Bet (y 


Disregarding signs, it follows from equation (10) that 


Aaa (1+ =). 


MN [yy 


For the first approximation we may take 


es. (12) 


j 


and for geometrical reasons y = 2.5/4, say that the second approximation 
would be 


Aah (14+ (13) 


Let us integrate differential equation (3a): 
U=—-2H,«+C, 


remembering U is magnetic potential difference between the pole cross-sections - 
symmetrically placed with respect to the medial plane, and not a magnetic 
potential. The constant of integration C is determined by the boundary con- 
dition U = U,, x = d,. We then take U, and H,, from (5) and (11) with approxi- 
mation (12) to obtain the magnetic potential difference distribution 


igre. eas 
Ue oe (1 N a1: (14) 
Substituting expressions (2), (14) and @ = 7° z B,, in the differential equation 
(1a), we find after some rearranging 


m E/e (1 + (6 — &)/A) 
V& + o2 arsinh (€/e) — m (1 + (0 — €)/A) 


We | 5D 


(15) 


This differential equation surely has no elementary solution; for the given 
parameters m, 0,4 and boundary condition &= 06; @=0, we must carry out 
numerical integration. The Runge-Kutta method, which is exceedingly precise 
for this calculation by virtue of the approximations made, yields with para- 
meters m = 1.25;6 =6;A= 21; Q: = 11 the polepiece profile curve plotted in 
Figure 7 as the solution of the differential equation (15). The slope of the 
tangent at the point = 6; 9 = 0, of a polepiece profile curve is obviously 
important: 


fn) 
ky —— = a oe 
(g:/m) V0? + of arsinh (/g,) — 0, 
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_ By use of the function x = x(n) defined by (7), it can be modified 


v mx o 

_At the point = 6 + A; @ = 9, the tangent is horizontal. 

Engineering practice will probably adopt some simplifications, rejecting 
curved tapering of the polepieces. In the present paper we shall assume a pole- 

_ piece shape already proposed by ANDREW and RusHworrtu [1], namely, conical 

tapering at the gap-end of the otherwise cylindrical polepieces. If the normal- 

ized radius g, of the cylindrical part is known, the frustum of the cone will be 
determined by the slope k, and the normalized radius 9, of its smaller base. Of 
course, it is not convenient to integrate the differential equation (15) for every 
particular case, to obtain radius @y. 

Taking into consideration the influence of the yoke upon the field distribu- 

tion, let us compare the fields depicted in Figures 1 and 8. Besides the visible 
difference, the field shown in Figure 8 is no longer rotationally symmetrical. 
Let us make a seemingly unacceptable assumption: that permeance of the 
elementary air shells does not depend on their position in space, but has values 
equal to that of the elementary air shell of the gap-end (more precisely, the 

initial air shell of the polepiece with neglected slope k,). Accordingly, assuming 
the same magnetic potential distribution as before and substituting 


oP 0/01 
dx Payiinde + (6/o,)2 arsinh (5/o,) 


= [oT % 


in the differential equation (1a) we obtain 


d 

Ee 0”) = Ws is 14 (1 +- 7 3 a . 

After integration that takes into account boundary condition € = 6, @ = Q, 
solution of this simple differential equation reads: 


e+ —* €-6-AP=olt min. 


Hence, polepieces are parts of the ellipsoids of revolution. Circular yoke-end of 
the polepiece with abscissa = 6 + A lies in the equatorial plane of the ellipsoid. 
Therefore 

oo = V+ mix. (17) 


From the above ‘exact’ formula we can derive the approximative one 


mh x 
= 1 
Ose Coches (18) 
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Referring to the numerical example, whose solution has been depicted in 
Figure 7, we are now able to estimate not only what we called the seemingly | 
unacceptable ‘equal permeance assumption’ concerning air shells, but also the” 
relative validity of the formula (18). The Runge-Kutta method and formulae 
(17) and (18) give the following values for the yoke-end normalized radius 

respectively: 11.98; 12.04; 12.1. We must conclude that air permeance distri- 
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bution is far less significant than magnetic potential difference distribution 


120 
— = 
iL 
: | 
| —- 
i= 
| 
Hot: 
| 
ee | a t le | 
110 == i 
6 10 14 18 22 26 
é ——— 
Figure 7 
Polepiece profile curve as the solution of the differential equation (15) for parameters given in 
the text. 


along the pole axis. This fact is the main justification for the usefulness of 
assuming circular lines of force (as was done in Section 2). Lastly by means of 
formulae (16) and (18) one can find the normalized height of the cone frustum 
(see Figure 1 and equations (6)): 


(1-7). (19) 


5. Optimal Design 


Given a set of initial conditions, e. g. induction in the gap and the gap 
volume, or induction in the gap and the gap ratio Qo, One can consider the 
question concerning the minimal volume V of the polepiece, once the magnetic 
material has been chosen. 


In general, we should look for the particular polepiece profile function 
y = r(x), which will make the volume 


ay td 
V=% / Yas 
ay 
as small as possible. This poses a very complic 


all as po ated problem in variation theory. 
Bearing in mind the lack of information concerning magnetic material, and also 
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for other reasons, we shall restate the optimal design problem thus: given a set 
of initial conditions, we must find the operating point on the demagnetization 
curve of the chosen magnetic material that will make the volume of the pole- 
‘piece as small as possible, provided that the design calculation follows the lines 
laid down in Section 4! 

Consider the problem with fixed B,, 7) and dy. This implies that @) is known. 
As before, we shall define the normalized polepiece volume 


V 


aesge 3 


“which can be calculated from the following somewhat cumbersome formula 


(m i x)8 


o 2 3A , 1 A 5 
It is interesting to note that 
B2 
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according to equations (4), (11) and (12). Hence, if the product B,, H,, is a 
maximum, the product m A will be a minimum. But as long as the coefficients 
of the polynomial (20) depend on the parameters 6, @,, x and A, the minimum 
of o does not necessarily coincide with minimum of mA. Taking the fixed 
operating point on the demagnetization curve, i.e. m and A as constant, we 
are now in a position to examine whether the minimum points of such curves 
as are depicted in Figure 5 coincide with the minimum of o. Since at these 


points ; 
dQ, _,. 4% 
oe ee 
we find that 
do 
75 ae 


By this account the minimal value of « becomes shifted towards a larger value 

of 6. Fortunately, numerical examples show that displacements of minimum 
points are entirely negligible. i 

Suppose the initial conditions B, and g, are given. Design procedure is then 
as follows: 

Let us choose an operating point on the demagnetization curve for the 
selected magnetic material. Thus parameters m and A are determined. Such 
curves as depicted in Figure 6 (if constructed for given @9) enable us to find 
optimal values for 6 and 9. The corresponding value of x can then be derived 
from Figure 3. However, at the first approximation, parameter x may be 
replaced by 1. Formula (20) now yields the normalized volume o. By the same 
means we can calculate values of o for several operating points on the demagnet- 
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ization curve in the neighbourhood of the (BH),,,, point, and plot the function 
o = o(m). Finally, minimum of this function determines the optimal induction 
ratio m. As is shown in the next Section, the optimal operating point will be 
found close to (BH),,,,, point for all practical cases. Having the optimal induc- 
tion ratio m, we can find the remaining parameters of the magnet. 


6. Verification and Application of the Developed Theory . 


Example 1. Permanent magnet installed at the University of St. Andrews, | 
Scotland, in 1951. ANDREW and RusHWoRTH [1] have given detailed informa- | 
tion about their permanent magnet, made from Alcomax III magnetic material. — 
Their initial conditions would have to be B, = 6000 G, 9, = 3.7. Reduced 
demagnetization curve for Alcomax III has been shown in Figure 2. At 
(BH) pax Point m = 0.766; A = 15. Normalized volumes of the polepieces in the — 
neighbourhood of this point were calculated and the result plotted in Figure 9. _ 
The optimal operating point is at m = 0.76 where o = 4782. The shift of the — 
induction ratio m being less than 0.8°%, the magnetic state of the material will _ 
remain at or very near (BH),,,, point. From Figure 6 we find 6 = 1.95: 0, = 4.8, 
and by means of formulae (18) and (19) 0, = 6; e = 4.65. Table 1 shows discre- 
pancies between calculated and planned poles. Zero discrepancy of 6 is inter- 
esting, whilst other discrepancies are justified mainly by reduction of demagnet- 
ization curve, which phenomenon the authors of [1] discovered after construc- 
tion of their magnet. 


Table 1 
0) Or A 05 € o 
calculated pole . . 1.95 4.8 15 6 4.65 478 x 
planned pole .. . 1.95 4.44 se7 SS }0) 4.57 395 x 
%-discrepancy with 
planned pole .. . 0 8.1 9.5 8.1 175 21 


To verify derived formulae, we must use quantities obtained by measure- 
ments on the made magnet. Initial conditions then read: B,=5500 G, o,=3.61 
Now we are not looking for the optimal polecap, but wish to approach as 
closely as possible to actual polecap dimensions. The following parameters at 
(BH) nay point can be regarded as independent: m = 0.702: 6 = 1.9305 =73.018 
and these parameters with equation (9) yield the polecap slope ky = 0.865, so 
that 0, = 4.39. Equation (13) with n = 1, My, = 1000 gives 2 = 13.88. Remai- 
ning parameters were calculated for x = 0,95: 0, = 5.44 and e = 4.51. Table 2 
shows discrepancies between calculated and made poles. 
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Figure 9 Figure 10 
Plot of the function o = o(m) giving the Plot of the function o = a(m) giving the 
optimal induction ratio m for the case of optimal induction ratio m for the case of 
Example 1. Example 2. 
Table 2 
Pee 
calculated pole. . . 4.39 13.88 5.44 4.51 
made Polen <7. «a. 4.33 13.36 5.42 4.45 


%-discrepancy with 
made-pole ... . 1.4 si) 0.4 1.4 


EE EEE a ee 
Except for the parameter ¢, positive discrepancies are due to oversizing of the 


poles, made throughout the design calculation. 
) 


Example 2. The proposed magnet. A large Columax magnet is to be designed 
with initial conditions B, = 11000 G, @ = 7. 
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The reduced demagnetization curve for Columax has been shown in Figure 2. 
At (BH),,qy point m = 1.25, A = 21.6. Normalized volumes of the polepieces in 
the neighbourhood of this point were calculated and the result has been plotted 
in Figure 10. The optimal operating point is at m = 1.225, where o = 30602. 
The shift of the induction ratio m being 2%, the (BH),,,, point is maintained. 
From Figure 6 we find 6 = 7.45; 9, = 11.07 and by means of formulae (18) and 
(19) @. = 12.14; e = 9.2. The University of Leeds Columax magnet has the gap 
length 2d, = 1.4", volume of the polepiece something about 41.1 dm3, and 
initial conditions B, = 9500 G, oy = 7.15. With our initial conditions and the 
above gap length, the proposed magnet would have a volume of about 
54.2 dm?, 32% larger than the University of Leeds magnet. 


7. Conclusion 


From the preceding two examples it is possible to derive an interesting 
plausible law: if the induction in the air gap is to be increased for x % with 
respect to an already constructed magnet, the volume of magnetic material would 
be increased roughly by 2 x %, provided that the gap ratio remains unaltered 
and that magnets are of optimal dimensions. 

As has been frequently emphasized, this method of calculation leads to 
oversizing of the magnet. This can be turned to good account, for it is well 
known that after magnetization permanent magnets are usually stabilized by 
means of a certain demagnetizing field. 

At all events, acceptability of the developed method will be established only 
after numerous applications. 
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Zusammenfassung 


In dieser Arbeit wird eine Methode zur Berechnung von grossen Permanent- 
magneten mit rotationssymmetrischen Polen entwickelt. Zum Unterschied von der 
klassischen Methode, die die Hopkinsonschen Streufaktoren benutzt, gibt diese 
Methode einen besseren Einblick in die geometrischen Verhdltnisse und fiihrt zu 
einfachen Formeln fiir die Berechnung der Magnetpoldimensionen. 
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Die wesentlichen Resultate sind: 

1. Bei gegebenen Luftspaltdimensionen und gegebener Luftspaltinduktion be- 
stehen bei vorgezeichnetem Magnetmaterial optimale Polhorndimensionen. 

2. Die Verteilung des magnetischen Skalarpotentials langs der Magnetpolachse 
ist weit wichtiger als die Magnetfeldverteilung in der Luft. 

3. Der giinstigste Betriebspunkt auf der Entmagnetisierungskurve ist gegeniiber 


dem Punkt (B H),,,, nur fiir wenige Prozente nach grésseren Werten von B ver- 
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Uber den Zusammenhang zwischen der Fliessbedingung 
eines elastisch-plastischen Korpers und seinem Fliessgesetz 
Von Hans ZIEGLER, ETH, Ziirich 


1. Einleitung 


Die von v. MisEs [1]4) aufgestellte Theorie des plastischen Potentials kann in 


_ der von PRAGER [2] verallgemeinerten Gestalt etwa wie folgt formuliert werden: 


Es seien g, (k = 1, 2, ..., 2) die verallgemeinerten Verformungen des ela- 
stisch-plastischen Kérpers, das heisst also ein Satz von Lagekoordinaten, 
welcher seine Deformation mit der gewiinschten Genauigkeit beschreibt. Ihre 


 elastischen und plastischen Anteile seien mit gj bzw. gf bezeichnet. Ferner sei 


(unter Verwendung der bekannten Summationsregel) 


dA = Q, day (1.1) 


die elementare Deformationsarbeit. Dann kann man die Q, (k = 1, 2,..., ”) 
als verallgemeinerte Krafte bzw. nach PRAGER als verallgemeinerte Spannungen 
bezeichnen. Sowohl der Verformungs- wie der Spannungszustand lassen sich 
durch Punkte mit den Fahrstrahlen q bzw. Q im n-dimensionalen euklidischen 
Spannungsraum R,, darstellen. Definiert man hier zum Zwecke der Metrisierung 


- das Skalarprodukt zweier Vektoren mit 


Q°q=Q2x%» (1.2) 
so ist die elementare Deformationsarbeit (1.1) durch 
dA=Q-:-dq (153) 


gegeben. Es entspricht dann jedem verallgemeinerten Verformungs- bzw. 


1) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 68, 
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Spannungszustand q, bzw. Q, ein Vektor q bzw. Q oder, wie wir im folgenden 
auch sagen werden, ein Punkt q bzw. Q in R, und umgekehrt. 

Betrachtet man den K6rper in einem gegebenen Stadium des Verformungs- 
prozesses, so lassen sich vom momentanen Spannungszustand aus gewisse 
andere ohne plastischen Fluss realisieren. Diese sollen als aplastische Spannungs- 
zustinde und die Gesamtheit ihrer Bildpunkte in R, als aplastischer Bereich R,, 
bezeichnet werden. Ferner sei unter einem aflastischen Spannungsinkrement 
dQ* ein solches verstanden, das zwei benachbarte aplastische Spannungszu- 
stande Q und Q’ ineinander iiberfiihrt. Die Fliessgrenze, eine Hyperflache in 
R,, wird dann als Oberflache des aplastischen Bereiches durch diejenigen apla- 
stischen Punkte Q gebildet, fiir die nicht jedes hinreichend kleine Spannungs- 
inkrement dQ aplastisch ist. Liegen die benachbarten Punkte Q und Q’ an 
der Fliessgrenze, dann nennt man das aplastische Spannungsinkrement dQ* 
neutral. Liegt von den beiden Punkten genau einer an der Fliessgrenze, so 
stellt dQ* eine Entlastung bzw. Belastung dar. 

Der Zusammenhang zwischen dem Spannungszustand Q und dem elasti- 
schen Verformungszustand q? wird wie in der Elastizitatstheorie angenommen, 
Fiir die im aplastischen Bereich R,, aber nicht an der Fliessgrenze liegenden 
Spannungspunkte ist das plastische Verformungsinkrement dq? null. Fiir die 
Punkte an der Fliessgrenze kann es gleich oder ungleich null sein. Beziiglich 
seiner Richtung gelten folgende Satze: 

Satz 1: Mit Q’ und Q” ist auch jeder Spannungszustand 


Q=0'+.2(Q"— 0), <aat) (1.4) 


aplastisch. Mit anderen Worten: Mit zwei beliebigen Punkten Q’ und Q” liegt 
auch ihre ganze Verbindungsstrecke im aplastischen Bereich, das heisst der 
aplastische Bereich ist konvex. 

Satz 2: Ist Q ein Spannungszustand an der F lessgrenze und dQ* ein belie- 
biges aplastisches Spannungsinkrement, das vom wirklichen Spannungsinkre- 
ment dQ véllig unabhangig ist, so gilt 


dQ*- dq’ <0. (1.5) 

Mit anderen Worten: Die Normalprojektion des plastischen Verformungs- 
inkrements auf jedes mit dem wirklichen Spannungszustand vertragliche apla- 
stische Spannungsinkrement ist nichtpositiv. 
Satz 3: Fiir das wirkliche Spannungsinkrement dQ gilt 


GQ agree t, (1.6) 


Das heisst : Die Normalprojektion des plastischen Verformungsinkrements auf 
das wirkliche Spannungsinkrement ist nichtnegatiy, 


2. Problemstellung 


| R. v. Mises hat den durch die Satze 1 bis 3 beschriebenen Sachverhalt 
oo als Hypothese und zudem unter Beschrankung auf ein Raum- 
element ausgesprochen, indem er als Verformungen q, und Spannungen Q,, die 
lokalen Verzerrungskomponenten ¢;; und die Spannungskomponenten o;; 
_ beniitzte. Es ist das Verdienst von PRAGER, die Theorie im Sinne von Abschnitt 
al verallgemeinert zu haben. Ihre Begriindung ist auf verschiedenen Wegen 
geurch BisHor und Hitt [2a], DRucKER [3] und den Verfasser [4] versucht 
worden, vermag aber insofern vielleicht nicht restlos zu befriedigen, als in allen 
_ Fallen gewisse physikalische Postulate beniitzt werden mussten. Der Verfasser 
: hat sich daher neuerdings die Aufgabe gestellt, unter Verzicht auf alle physi- 
_kalischen Annahmen wenigstens in aller Strenge nachzuweisen, dass aus der 
_ Voraussetzung, die Sitze 1 bis 3 seien in der v. Misesschen Fassung, das heisst 
J fiir das Raumelement richtig, auch ihre Giiltigkeit in der Pragerschen Verall- 
_ gemeinerung, das heisst im Sinne von Abschnitt 1 folgt. 
d Obschon die Theorie meist fiir elastisch-plastische Kérper formuliert zu 
werden pflegt, beschrankt man sich bei ihrer Anwendung gewohnlich auf den 
einfacheren starr-plastischen Kérper, der durch die Bedingung q’ = 0 bzw. 
q=¢ definiert wird. Fiir diesen K6rper ist der erwahnte Nachweis in einer 
kiirzlich verdffentlichten Arbeit [5], an die sich die folgenden Ausfiihrungen 
 anschliessen, erbracht worden. Hier handelt es sich darum, den Beweis auf 
- elastisch-plastische Korper zu tibertragen. Das ist unter der Voraussetzung 
_ méglich, dass die Aufspaltung des Verformungszustandes q in die Anteile q° 
und g? in geeigneter Weise vorgenommen wird. Freilich ist das Problem damit 
noch nicht so weit vereinfacht wie im Falle des starr-plastischen Kérpers, da 
sich im allgemeinen mit fortschreitendem plastischem Fluss der aplastische 
Bereich auch bei elastisch-idealplastischem Material verformt und damit ein 
4hnliches Verhalten zeigt wie bei der Verfestigung eines starr-plastischen 
_ Korpers. 

Die Verformungen werden im folgenden einfachheitshalber als klein und 
der Zusammenhang zwischen ihren elastischen Anteilen sowie den Spannungen 
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_als linear angenommen. 


3. Verformungen und Spannungen 


An einem Element des betrachteten Korpers sind der Verformungs- und der 
Spannungszustand durch die Verzerrungskomponenten ¢; ; bzw. die Spannungs- 
_komponenten o,;; gegeben; sie lassen sich durch je einen Vektor e bzw. s im 
neundimensionalen euklidischen Spannungsraum R, darstellen. Definiert man 


hier das Skalarprodukt mit 
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so ist die elementare Deformationsarbeit je Raumeinheit durch 


gegeben. Mit Riicksicht auf die Symmetrie der Tensoren ¢;; und o;, liegen die 
Vektoren e und s in Wirklichkeit in einem sechsdimensionalen linearen Unter- 
raum FR, von Ry. Hier entspricht jedem Verzerrungs- bzw. Spannungszustand 
&;; bzw. o;; ein Vektor e bzw. s und umgekehrt, und der aplastische Bereich 
wird wie in Abschnitt 1 definiert. Er kann sich ins Unendliche erstrecken; so 
schliesst man ja gewohnlich plastischen Fluss unter einem hydrostatischen 
Spannungszustand aus. 

Der Verformungs- und der Spannungszustand des ganzen Kérpers werden 
durch die von den Koordinaten x; abhangigen Funktionen e, ; bzw.o;, dargestellt, 
und diese lassen sich, wie SYNGE [6] und andere Autoren gezeigt haben, als 
Vektoren E bzw. S im Funktionenraum F interpretieren. Da die elementare 
Deformationsarbeit durch das tiber den ganzen Kérper K erstreckte Raum- 
integral 


th [o ; de,; dv (3.3) 
K 


gegeben ist, empfiehlt es sich, fiir die Metrisierung des Funktionenraums das 
Skalarprodukt durch 


S-E=/o,;6,;dv (3.4) 
x 


zu definieren. Wie in [5] gezeigt wurde, ist diese Definition zulassig; sie ergibt 
eine positiv definite Metrik sowie fiir die elementare Deformationsarbeit (3. 3) 
die Darstellung 


dA =S-dE. (3.5) 


Jedem Verformungs- bzw. Spannungszustand im K6rper K entspricht ein 
Vektor im Funktionenraum. Dabei beschrinkt sich diese Darstellung keineswegs 
auf kinematisch zulassige Verformungs- und statisch zuldssige Spannungszu- 
stande; vielmehr sollen vorderhand alle denkbaren Zustande zugelassen sein, 
fiir welche die Integrale (3.3) bzw. (3.4) existieren, Der aplastische Bereich F 
kann sich auch hier ins Unendliche erstrecken. Lasst man die Belastung des 
Kérpers K von Null an wachsen, dann durchwandert der Bildpunkt § des 
wirklichen Spannungszustandes zunachst den aplastischen Bereich und 
erreicht die Fliessgrenze im Augenblick, in dem die ersten Elemente von K 
plastisch werden. Mit weiter zunehmender Belastung setzt in einer sich all- 
mahlich erweiternden Umgebung dieser Elemente plastischer Fluss ein, und 
es kann gegebenenfalls zum Kollaps kommen, ohne dass aber der ganze K6érper 
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ye 


_ plastisch zu werden braucht. Ein wesentlicher Unterschied zum starr-plasti- 
_ schen KG6rper besteht darin, dass hier der plastische Fluss lokal einsetzen kann, 
-sobald die ersten Elemente an ihrer Fliessgrenze angelangt sind, wahrend dort 
die starre Umgebung des plastifizierten Gebietes den lokalen Fluss verhindert 
und plastische Verformungen erst auftreten, wenn ein geniigend grosser Teil- 
bereich von K plastifiziert ist. Hieraus folgt fiir die Darstellung im Funktionen- 
raum 
Satz 4: Der aplastische Bereich des ganzen elastisch-plastischen Kérpers 

deckt sich im allgemeinen nicht mit demjenigen des entsprechenden, lokal der 
gleichen Fliessbedingung unterworfenen starr-plastischen K6rpers. 
; Insbesondere sieht man ohne weiteres ein, dass der aplastische Bereich des 
_ jungfraulichen elastisch-plastischen Kérpers in demjenigen des entsprechenden 
_ starr-plastischen Kérpers enthalten ist. Aus der Méglichkeit des lokalen 
_ Flusses ergibt sich beilaufig auch eine Aussage fiir den idealplastischen, das 
- heisst den Korper ohne Verfestigung, dass naémlich im elastisch-idealplastischen 
é im Unterschied zum starr-idealplastischen Kérper Spannungszustande reali- 
 sierbar sind, deren Bildpunkte ausserhalb der urspriinglichen Fliessgrenze 
- liegen. Somit hat man 
Satz 5: Im Gegensatz zum aplastischen Bereich des starr-idealplastischen 
_ Korpers ist derjenige des elastisch-idealplastischen K6rpers nicht fest, sondern 
_ andert sich im allgemeinen, wahrend die Belastung bis zum Kollaps gesteigert 
_ wird. 
: Der elastisch-idealplastische Koérper verhalt sich also mit zunehmender 

Ausdehnung des plastifizierten Gebietes ahnlich wie der starr-plastische Kérper 
- bei der Verfestigung. 

Der Ubergang zu verallgemeinerten Verformungen und Spannungen wird 

_ stets dadurch vollzogen, dass man sich auf die Betrachtung von speziellen Ver- 
_ formungszustanden beschrankt, die durch eine endliche Zahl n von Parametern 

g, beschrieben werden kénnen. Man kann sich diese Restriktion in der Weise 

durchgefiihrt denken, dass man geeignete Bindungen einfiihrt (und allenfalls 
- auch vorhandene Bindungen lést), derart, dass sich der Freiheitsgrad des 

Kérpers (von einer allfalligen starren Bewegung abgesehen) auf  reduziert. 
Es ist klar, dass dieser Prozess einen groben Eingriff darstellt, der nur dann 
zu verantworten ist, wenn er das wirkliche Verhalten des Korpers in der gege- 
 benen Situation nicht zu stark modifiziert, der aber unter dieser Voraussetzung 
das Problem wesentlich vereinfacht. 
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4. Beispiel 


. Es sei ein Balken mit rechteckigem Querschnitt unter Beanspruchung auf 
_ Zug sowie spezielle Biegung betrachtet und hievon in Figur 1 eine Scheibe 
der Dicke 1 ins Auge gefasst. Der Verformungs- und der Spannungszustand 
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dieser Scheibe werden, falls man die Langeneinheit hinreichend klein gewahlt 
hat, durch je sechs von y und z abhangige Verzerrungs- bzw. Spannungskom- 
ponenten beschrieben. Man kann aber das Problem im Sinne der Festigkeits- 
lehre durch die Annahmen vereinfachen, dass die Querschnitte bei der Defor- 
mation in erster Naherung eben bleiben und dass sich die Langsfasern nicht 
wesentlich beeinflussen. Das heisst, man kann sich vorstellen, dass der Zu- 
sammenhang zwischen den Langsfasern geloést (der Spannungszustand mithin 
linear) und die Bewegung jedes Querschnittes diejenige einer starren Flache 
sei. Wie Figur 2 zeigt, wird die Verformung der Scheibe in dieser Naherung 


i 


4 


Figur 1 


Figur 2 
Balkenelement unter Beanspruchung auf Zug 
und Biegung. 


Verformung des Elements. 
durch die Dehnung ¢ sowie die Kriimmung der elastischen Linie, nimlich den 
Winkel x zwischen den beiden Stirnflachen beschrieben, und dabei gilt 


e=7%. (4.1) 


Die zugehorigen verallgemeinerten Spannungen sind die Normalkraft N sowie 
das Biegemoment MM, da die elementare Deformationsarbeit durch 


dA =Nde+M dx 
gegeben ist. 

Um den urspriinglichen aplastischen Bereich der von jetzt an als elastisch- 
idealplastisch angenommenen Scheibe zu erhalten, hat man die Spannungszu- 
_ stande zu diskutieren, die an der Fliessgrenze liegen. Fiir den in Figur 3 dar- 
gestellten Fall, bei dem die Fliessgrenze O» in den obersten Fasern auf Zug 
erreicht wird, ist die Spannungsverteilung durch 


i alee 4 
oF ERE a = 
gegeben, die Normalkraft also durch 
: h/2 
eth eer 


—h/2 
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und das Biegemoment durch 


h/2 
Matlacdi=— 4.4 
as +h 70° (4.4) 
—h/2 


Wenn man aus (4.3) sowie (4.4) 7 eliminiert und die Falle hinzunimmt, in 
denen die Fliessgrenze auf Druck bzw. am unteren Querschnittsrand erreicht 
wird, erhalt man die Grenzen 


M=+4 4 (ho +N) (4.5) 


hs 


eines Parallelogramms, das in Figur 4 mit 0 bezeichnet ist und innerhalb des 
von PRAGER (vgl. [2], S. 46) gegebenen, durch zwei Parabeln begrenzten apla- 
stischen Bereiches fiir starr-idealplastisches Material liegt. 
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Figur 3 Figur 4 

; Spannungsverteilung bei Erreichen der Sukzessive Fliessgrenzen des plastisch 

o4 

i Fliessgrenze am oberenRand. gebogenen Elements. 


, Beansprucht man jetzt die Scheibe in reiner Biegung durch ein positives 
- Moment M, iiber die Fliessgrenze hinaus, so ergibt sich durch Anwendung der 
- Integrale (4.3) sowie (4.4) auf die Spannungsverteilung von Figur 5 der Zu- 


_ sammenhang 


M, = = (3h® 42) 6 (4.6) 


- zwischen der halben Hohe ¢ des elastischen Kerns und dem Biegemoment M,. 
- Bei elastischem Verhalten ware das Biegemoment fiir gleiche Deformation > 


D bh 
M, ieee (4. 7) 
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der plastische Fluss in den ausseren Fasern bedeutet einen Spannungsabfall 
gegentiber der elastischen Lésung, und diesem entspricht eine Reduktion des 
Biegemomentes um 


AM =M, = M, == ia ae ee (4.8) 


Figur 5 


Spannungsverteilung unter plastischer Biegung. 


Um den neuen aplastischen Bereich der durch M, auf Biegung vorbean- 
spruchten und dadurch einer ersten plastischen Verformung unterzogenen 
Scheibe zu erhalten, hat man jetzt Spannungszustande zu betrachten, die 
sich fiir |z|>|¢| durch den gleichen marginalen Spannungsabfall von der 
linearen unterscheiden wie in Figur 5. Unter diesen sind wiederum diejeni- 
gen zu diskutieren, die wie die Spannungszustande mit o(¢) = 0, bzw. 
o(— h/2) = og von Figur 6 und 7 an der Fliessgrenze liegen. Dabei kann man 


Figur 6 Figur 7 
Spannungsverteilung des plastisch gebogenen Spannungsverteilung des plastisch gebogenen 
Elements fiir o(f) = op. Elements fiir o( — h/2) = op. 


die Integrale (4.3) und (4.4) auf die elastischen Verteilungen anwenden und 
hat nur nachtraglich das Biegemoment um den Betrag (4.8) zu reduzieren. Im 
Fall von Figur 6 ergibt sich so, wenn man von der Spannungsverteilung 
y+ z 
Co = Wat ory (4.9) 
ausgeht, 


—- 69 -AM , (4.10) 
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a 

im Fall von Figur 7 mit 

E 28 =i 

dj Co = Ct (ry — hj2) (7 + 2) 09 (4.11) 
analog 

: 

ai een Ny: bh 2t—h/2 

© N=bhr——__ = 

j Flr — hj2) 7 sowie M 12 Er —h/2) % AM. (4.12) 

: 

A ‘Durch Elimination von 7 aus (4.10) bzw. (4.12), Hinzunahme der Falle 


(— ¢) =— 9, o(h/2) = — oy und Beriicksichtigung von (4.8) erhalt man 
_schliesslich die Grenzen 


h? 
- sowie 
i b h? 
M soare Lect Ges M,- a oP (4.14) 


_ des neuen aplastischen Bereiches, der wiederum die Form eines (beziiglich der 
gz Achse M symmetrischen) Vierecks besitzt. Er ist in Figur 4 fiir € = h/4 bzw. 
_€ =/h/10 eingetragen und mit 1 bzw. 2 bezeichnet; ferner ist leicht einzusehen, 
_ dass er sich fiir € > 0 auf die Teilstrecke 3 


; vss, (4.15) 
der Achse M zusammenzieht. 

Wie zu erwarten war, folgt aus diesem Beispiel 

Satz 6: Die Satze 4 und 5 gelten nicht nur fiir die wirklichen, sondern auch 

fiir die verallgemeinerten Spannungen des endlichen K6rpers. 

Figur 4 zeigt deutlich den Unterschied, der zwischen den aplastischen 

_Bereichen des elastisch-plastischen und des starr-plastischen K6rpers besteht. 

_Freilich behalt wenigstens fiir den elastisch-idealplastischen K6érper der apla- 

stische Bereich des starr-idealplastischen Kérpers seine Bedeutung insofern, 

als auch hier mit dem Kollaps zu rechnen ist, sobald der Spannungspunkt seine 

_Grenze erreicht. 


5. Elastische und plastische Verformungen 


Fiir die Theorie des plastischen Potentials ist die Art der Aufspaltung des 
Verformungszustandes in seinen elastischen und den plastischen Anteil von 
grundlegender Bedeutung. Hier liegt ein Problem, das in der Literatur meist 
iibergangen, tatsdchlich aber nicht ganz trivial ist. 
_ Beim Raumelement pflegt man, vom linearen Spannungszustand mit dem 
_Spannungs-Dehnungs-Diagramm von Abbildung 8 ausgehend, denjenigen An- 
teil ef, des Verzerrungszustandes ¢;; als elastisch zu bezeichnen, der dem Span- 
| Weinecrustand o,; unter der Annahme rein elastischen Verhaltens, das heisst 
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nach dem Hookeschen Gesetz entspricht; der plastische Verzerrungszustand 
ist dann durch die Differenz ef, = e;; — &; definiert. Im urspriinglichen apla- 
" stischen Bereich gilt «?, = 0 baw. e;; = e&{;; ferner sind die ef, diejenigen Ver- 
zerrungen, die sich bei der Entlastung zuriickbilden, wahrend die e?, bestehen 
bleiben, und schliesslich stellt im Ausdruck (3.2) fiir die elementare Deforma- 
tionsarbeit je Raumeinheit, der jetzt mit 


angeschrieben werden kann, der erste Term rechterhand die Zunahme der 
elastischen Verzerrungsenergie und der zweite die elementare, in Warme tiber- 
gefiihrte Dissipationsarbeit dar. 

Es sei nur nebenbei vermerkt, dass die zuletzt genannten Eigenschaften 
nicht ohne weiteres aus der vorstehenden Definition der Aufspaltung folgen. 
So muss man, um ihre Aquivalenz auch bei grésseren Verzerrungen zu erhalten, 
fiir die Entlastung etwa vom Punkt C in Figur 8 den Weg des Bildpunktes mit 
CDE und nicht mit CDF, den Bauschingereffekt mithin als beschrankt an- 
nehmen. Um die Schwierigkeiten zu vermeiden, die sich alsdann bei einer Bela- 
stungsumkehr ergeben, berufen wir uns auf die bereits am Schluss von Abschnitt 
2 eingefiihrte Voraussetzung kleiner Deformationen. 


Figur 8 


Spannungs-Dehnungs-Diagramm eines elastisch-plastischen Materials mit Verfestigung. 


Fir den ganzen Kérper geht man zweckmassig vom Spannungszustand S 
aus und spaltet unter Beniitzung der zugehérigen Spannungskomponenten o; , 
in jedem Punkt den Verzerrungszustand €;; m der eben besprochenen Weise 
auf. Damit hat man im urspriinglichen aplastischen Bereich wieder E’ = 0 bzw. 
E = E*. Nach dem Einsetzen der plastischen Verformung stimmt aber der 
Spannungszustand S im allgemeinen nicht mehr mit demjenigen §* iiberein 
der in K fiir die gegebene Lagerung und Belastung bei rein elastischem Ver 
halten zu erwarten wire. Bei der Entlastung bildet sich § im Gegensatz zu S* 
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im allgemeinen nicht zuriick; die Verformungen E* und E? lassen sich also 

> als diejenigen Anteile des Verformungszustandes deuten, die bei der 

- Entlastung von  verschwinden bzw. zuriickbleiben. Dafiir folgt aber bei 

a Art der Aufspaltung aus (3.3) und (3.5) die elementare Deformations- 
-arbeit 

dA =S-dE°+S- dE?, (52) 


und zwar stellt nach (5.1) der erste Term rechterhand die Zunahme der ela- 
stischen Deformationsenergie und der zweite die elementare Dissipationsarbeit 
dar. 

Wie am Schluss von Abschnitt 3 ausgefiihrt wurde, wird der Ubergang auf 
verallgemeinerte Verformungen und Spannungen durch Beschrankung auf Ver- 
formungszustande vollzogen, die sich durch Uberlagerung aus  speziellen Zu- 
standen gewinnen lassen. Im Funktionenraum bedeutet dies (vgl. [5]) stets die 
Restriktion auf einen linearen Unterraum F,, der den Ursprung enthalt und 
durch eine orthonormale Basis I, (k =1, 2, ... , m) dargestellt werden kann, so 
dass 


LS EE Cee EET 


oN S> 


E=iE- 1) I, (5.3) 


Noe Nee eon 


: ist. Die zugehérigen Spannungszusténde werden im allgemeinen nicht durch 
_ Vektoren S aus F, dargestellt. Indessen kann ein beliebiger Spannungszustand 
.S gemass 

J SaaS +-5,, S=(S- hy - (5.4) 

in seine Normalprojektion S$’ auf F, und einen zu F, orthogonalen Vektor S, 

_zerlegt werden. Nach (3.5), (5.3) und (5.4) ist dann die elementare Deforma- 

_ tionsarbeit durch 

c dA=S-dE=(S-1I,) dE-I,) =S’:-dE (ro) 


_ gegeben. Auf der Moglichkeit, diese Arbeit als Skalarprodukt zweier Vektoren 
‘im linearen Unterraum F, darzustellen, beruht beim starr-plastischen K6rper 
(E = E?) die Ubertragung der Satze 1 bis 3 vom Element auf verallgemeinerte 
_Verformungen und Spannungen. 

Fiir die Aufspaltung des Verformungszustandes kann man auch hier — dies- 
mal unter Beschrankung auf Verformungszustande E in F, — vom Spannungs- 
-zustand § ausgehen und unter Beniitzung der Spannungskomponenten o;, den 

- Verzerrungszustand ¢,; in jedem Punkt in seinen elastischen bzw. plastischen 
_ Anteil auflésen. Da aie mit dem Verformungszustand E seine Anteile E* und 
_ E? nicht notwendig im linearen Unterraum F, liegen, ist die Behandlung in 
-verallgemeinerten Gréssen auf dieser Basis im allgemeinen nicht médglich. 
Vielmehr empfiehlt sich, von den verallgemeinerten Verformungen gq und den 
-zugehérigen verallgemeinerten Spannungen Q auszugehen und die Zerlegung 


, 
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auf Grund der Forderung vorzunehmen, dass in dem aus (1.3) folgenden Ausdruck 
dA =Q-dq’ + Q dq? (5.6) 


fiir die elementare Deformationsarbeit der erste Term rechterhand die Zunahme 
der elastischen Deformationsenergie und der zweite die elementare Dissipations- 
arbeit darstellen soll?), Wie der Vergleich mit (5.2) zeigt, ist dann 


Q-dqe=S- dE =dU, (5.7) 


wobei U die elastische Deformationsenergie darstellt, sowie 


Q - dq’ =S - dE?, (5.8) 


und q? kann als derjenige verallgemeinerte Verformungszustand interpretiert 
werden, der dem verallgemeinerten Spannungszustand Q bei rein elastischem 
Verhalten entsprechen wiirde. 


6. Beispiel 


Bei der in Figur 1 wiedergegebenen, auf Zug und Biegung beanspruchten 
Scheibe ist die Behandlung der Querschnitte als starre Flachen mit der Be- 
schrankung auf Verformungszustande in einem zweidimensionalen Unterraum 
Fy, von F gleichbedeutend. Diese sind insbesondere durch eine lineare Vertei- 
lung der axialen Dehnung «, iiber den Querschnitt charakterisiert und lassen. 
sich durch die beiden Parameter ¢ sowie x darstellen. Die zugehérigen Span- 
nungszustande (vgl. Figur 5 bis 7) weisen im allgemeinen diesen linearen 
Verlauf nicht auf und gehéren somit dem Unterraum F, nicht an. 

Geht man von einer dieser Spannungsverteilungen o,(z) aus, und ermittelt 
man mit Hilfe des Hookeschen Gesetzes fiir jeden Wert von z die elastische 
Dehnung ¢%(z), so zeigt sich ohne weiteres, dass diese sowie die plastische Deh- 
nung ef(z) mit der totalen e,(z) im allgemeinen nicht linear in z ist. Demnach ist 
die Behandlung des Problems im linearen Unterraum F, auf dieser Basis aus- 
geschlossen. Dagegen kann man die elastischen Anteile der verallgemeinerten 
Verformungen durch die der Elastizitatstheorie entnommenen Beziehungen 


ait 12M 
a, a a 
: Ones aa Eb} (6.1) 


definieren und erhalt hieraus die plastischen Anteile ef = ¢ — &%, x? = x — x, 
Die lokale elastische Dehnung ist dann 


Fees 12M 
= ar tae (6. 2) 


*) Das ist die vom Verfasser bereits in [4] verwendete Definition. 


Vol. XII, 1961 Zusammenhang zwischen Fliessbedingung eines Kérpers und Fliessgesetz 65 


und die elastische Deformationsenergie berechnet sich hieraus erwartungs- 
~ gemiiss zu 
h/2 3 
N 12 M2 
ee 2 
x >| Lele) Nt 2 eae FE 2) 


aa 


a a 


= 


7. Der ganze K6rper 


Es sei jetzt angenommen, dass Satz 1 zwar fiir jedes Element des Korpers K, 
aber nicht fiir den ganzen Spannungszustand in K gelte. Dann gibt es zwei 
_ aplastische Spannungszustaénde §’ und S” sowie einen Wert a so, dass der 
_ Spannungszustand 


ee RN) OL 


Ses fais = 8S) (0a 1) (7.1) 


nicht aplastisch ist. Das ist aber nur méglich, wenn in wenigstens einem 
_ Element von K die Fliessgrenze iiberschritten, also hier im Widerspruch mit 
der eben getroffenen Voraussetzung mit o;; und of der Spannungszustand 
0;;=0;,+4(6,—0;), 0<a<1) (7x2) 
a 

4 nicht aplastisch ist. 
| Ebenso einfach ist der Nachweis, dass fiir einen Spannungszustand S$ an 
_ der Fliessgrenze, ein beliebiges aplastisches Spannungsinkrement dS* und das 
- wirkliche plastische Verformungsinkrement dE? die Ungleichung 


i dS* - dE? <0 (7.3) 
_ bzw. nach (3.4) | 
4 iy do*, de?, dv <0 (7.4) 


 besteht, falls fiir jedes Element von K Satz 2 richtig ist. Man kann naémlich 
~ diese Elemente in zwei Gruppen einteilen, je nachdem sie unter dem wirklichen 
_ Spannungsinkrement dS plastisch fliessen oder nicht. Fiir die Elemente der 
_zweiten Gruppe ist de?. = 0; sie liefern also keine Beitrage zum Integral (7.4). 


Fiir diejenigen der ersten Gruppe gilt nach Satz 2 
do;; de = eai0? (7.5) 


“ihre Beitrage zu (7.4) sind demnach nichtpositiv. 
Schliesslich folgt aus der Giiltigkeit von Satz 3 fiir jedes Element von K, 
dass fiir das wirkliche Spannungsinkrement dS die Ungleichung 


dS - dE’ =0 (7.6) 
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bzw. 
i do,, de? dv = 0 (7.7) 
K 


besteht. Die unter dS nicht fliessenden Elemente leisten naémlich keinen, die 
iibrigen nach Satz 3 nichtnegative Beitrage 

do,,dee 20 | (7.8) 
zum Integral (7.7). 


Somit hat man 

Satz 7: Gelten die Satze 1 bis 3 fiir das Raumelement, dann sind sie auch 
fiir die Darstellung des ganzen Verformungs- und Spannungszustandes im 
Funktionenraum richtig. 


8. VerallSemeinerte Gréssen 


Beim Ubergang auf verallgemeinerte Verformungen und Spannungen ist 
zu beachten, dass sich ein bestimmter verallgemeinerter Spannungszustand Q 
durch unendlich viele wirkliche Spannungszustande in K realisieren lasst, so 
dass der Punkt Q in R, unendlich vielen Punkten S in F entsprechen kann. 

In einem beliebigen Stadium des Verformungsprozesses fiihrt das verallge- 
meinerte aplastische Spannungsinkrement dQ* nach Abschnitt 1 zwei aplasti- 
sche Spannungszustande ineinander iiber und ist daher in R, Abbild eines 
aplastischen Spannungsinkrements dS* in F. Ferner liegt mit Q auch der wirk- 
liche Spannungszustand § an der Fliessgrenze, und nach (5.8) sowie (7.3) gilt, 
wenn dq’ bzw. dE? das wirkliche plastische Verformungsinkrement darstellt, 


dQ* - dq? = dS*- dE’ <0, (8.1) 


das heisst (1.5). Analog ergibt sich fiir das wirkliche Spannungsinkrement dQ 
bzw. dS * 

dQ -dq’=dS-dE’=>0, (8.2) 
somit (1.6). 

Sind ferner Q’ und Q" in R,, zwei verallgemeinerte aplastische Spannungs- 
zustande und S” sowie S” zwei ihnen entsprechende aplastische Zustande in F, 
so folgt daraus, dass die Spannungszustande (7.1) aplastisch sind, auch die 
Aplastizitat der entsprechenden verallgemeinerten Spannungszustinde Q. 
Stellen dq in R, und dE in F, dasselbe verallgemeinerte Verformungsinkrement 


dar, dann gilt nach (1.3) und (3.5) 
dA =Q-dq=S-dE, (8.3) 


und gemass Definition der verallgemeinerten Spannungen ergibt sich umge- 


a 
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kehrt aus dem Bestehen einer fiir jedes Paar dq, dE giiltigen Identitat von 
dieser Form, dass Q der zu S gehérende verallgemeinerte Spannungszustand 
ist. Nun folgt aus (8.3) fiir beliebige Paare dq, dE 


Q'-dq=S'-dE sowie Q"-dq=S"-dE (8.4) 


ee Pet eee 


und damit auch 


. 


[Q’ + 4 (Q" — Q’)] -dq = [S' + a(S" —S'))-dE. (8.5) 


Den Spannungszustanden S aus (7.1) entsprechen mithin die verallgemeinerten 
_ Spannungszustande 


CBU FAO — 0). 10 a= 1) (8.6) 


QP eee Va 


und daher ist mit Q’ und Q” auch jeder verallgemeinerte Spannungszustand 
_ (1.4) aplastisch. 


A Demnach hat man jetzt auch 

; Satz 8: Gelten die Satze 1 bis 3 fiir das Raumelement, dann sind sie auch 
_ fiir die Behandlung des ganzen K6rpers in verallgemeinerten Verformungen und 
_ Verzerrungen richtig. ; 


9. Bemerkungen 


Wie ohne weiteres aus ihrer Herleitung hervorgeht, beruhen die Satze 7 
und 8 auf der Voraussetzung, dass die Zerlegung des Verformungszustandes in 
seinen elastischen und plastischen Anteil nach den in Abschnitt 5 angegebenen 

_ Regeln vorgenommen worden ist. Es ist nicht schwer, die Giiltigkeit insbe- 
 sondere des letzten Satzes am Beispiel von Abschnitt 4 bzw. 6 zu tiberpriifen’). 
Die Tatsache, dass die Satze 1 bis 3 nicht nur fiir den starr-plastischen, son- 
-dern auch dem elastisch-plastischen Kérper gelten, darf indessen. nicht zur 
Annahme verleiten, dass die praktische Behandlung konkreter Probleme in den 
beiden Fallen gleich einfach sei. Die beiden Satze verkniipfen lediglich das ver- 
allgemeinerte plastische Verformungsinkrement mit der Fliessgrenze; diese 
selbst muss aber vermittelst einer Diskussion des ganzen Spannungszustandes 
S gefunden werden und lasst sich auf Grund der verallgemeinerten Spannungen 
Q allein nicht erhalten. Beim starr-idealplastischen K6rper ist das, wie ja schon 
das Beispiel in Abschnitt 4 zeigt, eine Aufgabe, die nur einmal zu lésen ist, da 
es hier nur einen einzigen, festen aplastischen Bereich gibt. Beim elastisch- 
idealplastischen Kérper dagegen deformiert sich der aplastische Bereich, und 


, 


3) Die hier unterdriickte Kontrolle ist mit Unterstiitzung durch die Herren Dipl. Math. W. Iss- 
LER und Dipl. Ing. H. M. FiscHER ausgefiihrt worden, denen der Verfasser auch die kritische 
Durchsicht des Manuskripts verdankt. 
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es muss daher nicht nur Schritt fiir Schritt das verallgemeinerte plastische 
Verformungsinkrement auf Grund der Satze 1 bis 3 neu ermittelt, sondern 
laufend unter Riickkehr zum wirklichen Spannungszustand auch die sich wan- 
delnde Form der Fliessgrenze erschlossen werden. Fiir den elastisch-plastischen 
Korper bedeutet daher die Behandlung in verallgemeinerten Gréssen lange 
nicht dieselbe Vereinfachung wie im starr-plastischen Fall. 


Abstract 


The author has recently shown that v. MisgEs’ theory of the plastic potential, if 
valid for the elements of a rigid-plastic body, also holds for the representation in 
generalized strains and stresses. This paper supplies the analogous proof for an 
elastic-plastic body. Here, however, the yield surface deforms, even in the case of 
perfect plasticity. Besides, it is impossible to determine the shape of the yield sur- 


face without considering the real state of stress. Thus, the simplification resulting — 
from the use of generalized coordinates is less than in the case of a rigid-plastic body. — 
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A Further Note on a Two-Dimensional Free-Boundary Problem 


By R.P.GILBErtT, East Lansing, Michigan, U.S.A.1) and 
H. KEnpDALL, Pittsburgh, Pennsylvania, U.S.A. ?) 


In a recent paper’), we obtained the following expression for the free-boundary 
between two fluids under dynamic conditions in a bent channel: 


y 
1 f ae ayylVeoa 
= 2, in Jeosh = Ks - zs = fe )| dy 
a VK? — K3 Cre! Ol =) | (1) 


Ve 
K 
_ (tanh5*) (y a Vo) ’ 


where 
K,=cosh[xcot(m— 9)], K,=sinh[mcot(m— g)], 
: 2; 
K, = cos (7 sin (x — 0) jena ii = eosh (7,”) ¢ 2) 
Voo Voo 


_ The boundary condition y = y, at ¥ = x, was specified by setting dv/dy = cot (p— It) 


at the point C; this gave 
J o-= 


Voo a ts (fie) 
uy Ki — Ki ad 
with %, y, cot(y — 2). 

This note is meant to serve as a preliminary report on the numerical evaluation 
of Equation (1), for various physical conditions. Since all of the computations in 
this note were performed on a desk calculator it was necessary to limit the range of 
physical parameters being used. More extensive computations are presently being 
made on an I.B.M. 702, and these will be reported on when completed. The 


following magnitudes were decided upon as being reasonable values for the physical 


and geometric parameters: 


p =2/2, the angle of the channel bend; 

f= 0°25, porosity; 

k = 100 milli darcys, permeability ; 

fy = 1 centipoise, viscosity of water; 

Voo = 15 ft/year 1-44 x 10-*cms-}, velocity of the water at oo; 
0» = 1:00 gcm-, density of water; 

0) = 0:90 gcm-, density of oil; 

y = (kg/tHw f) (@w — Qo) = 3°87 x 10-> cms”. 


A new parameter 6 = y 2/Voo = 8°5, and the dimensionless lengths % = «/a, y = y/a 
(where a is the channel width) were introduced because of their usefulness in the 


computations. 


I 


I 


1) Department of Mathematics, Michigan State University. 
2) Department of Mathematics, University of Pittsburgh. 
3) R. J. Bean, R. P. GILBERT, and H. Kenpat1, On a Two-Dimensional Free-Boundary Pro- 


blem, Z. angew. Math. Phys. 11, 341 (1960). 
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i i ; le employing 
Equation (1) has been evaluated by a two point SIMPSON s Tu nployir 
Reianie ad sizes; a check has been obtained by rewritting this equation in its 
derivative form and proceeding from point to point with an average slope. Below is 

a figure of both these curves. 


ih Ie 70 0-6 0-6 04 —*2 0 


Flow of water = sx 


= Sedat eS 0-2 
= Static oil-water 7 
infertace Saiulion using y 
Water SS Simpsons rule 0-4 
Solution using point ; 
ae ae S 0-6 
p=n/2 average slope method 
6 . 
08 
%-1-027-A7-9 
Zusammenfassung 


Diese Notiz enthalt eine numerische Illustration zur Verdffentlichung On a Two- 


Dimensional Free-Boundary Problem. Dabei wird vermittelst numerischer Integra- 
tion der in der erwahnten Schrift vorkommenden Differentialgleichung der Fall 
eines Kanals mit einer Kriimmung von 90° behandelt. Es ist beabsichtigt, zu einem 
spateren Zeitpunkt eine umfassendere Darstellung der Lésungen dieser Differential- 
gleichung zu veréffentlichen. 


(Received: September 12, 1960.) 


Ségrégation d’impuretés dans la glace 


Par CLAUDE JaccarD et LAURA Levi, Zurich?) 


1. Introduction et théorie 


Dans les travaux antérieurs [1, 2]*) nous avons examiné les propriétés électriques 
de la glace dotée d’acide fluorhydrique, et indiqué une méthode de production de 
specimens unicristallins. Toutefois lignorance du coefficient de ségrégation et des 
conditions précises de la croissance ne permettait pas d’obtenir des cristaux de 


concentration prescrite. C’est pour combler cette lacune que les recherches exposées 
ci-dessous ont été entreprises. 


1) Institut de Physique de PE.P.F., Zurich; 
neige et des avalanches, Weissfluhjoch/D 
Buenos Aires, Argentine. 

*) Les chiffres entre crochets renvoient al 


actuellement Institut fédéral pour l’étude de la 
avos, respectivement Instituto de Fisica de la Atmésfera, 


a Bibliographie, page.76. 
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1.1. Le probléme théorique est le méme que celui qui se pose pour la purification 
d’éléments semi-conducteurs et a déja été traité sous certains aspects dans des 
Ouvrages specialisés [3]. Cependant, comme la diffusion des impuretés est beaucoup 
’ plus faible dans l’eau a 0°C que dans les composés a point de fusion élevé, les con- 
‘ditions de la ségrégation effective sont différentes. Une solution approchée.a déja été 
‘calculée par Pout [4], mais elle n’est valable que pour une concentration initiale 
homogéne dans le liquide, et seulement une fois que la cristallisation a progressé 
d'une certaine distance; nous avons repris l'étude détaillée du phénoméne (dans un 
cas 4 une dimension) cherchant la concentration a n’importe quel endroit du cristal. 
En supposant que la solution s’étende jusqu’a l’infini, que la convection soit négli- 
geable et que le cristal croisse le long de l’abscisse x a une vitesse v constante, on 
cherche la répartition exacte des impuretés dans le liquide n(x, t) en connaissant le 
coefficient de diffusion D (pratiquement nul dans le solide) et le coefficient de 
ségrégation S = 1/s (rapport de la concentration de la solution v a celle de la glace 
a l’interphase n,). Pour cela il faut résoudre l’équation différentielle 


ESS 


avec la répartition initiale n(x, t = 0) = No + (x) et, en assignant l’abscisse vt a 
l'interphase, puisque les impuretés rejetées par le cristal sont emportées par le 
courant de diffusion: 


a ae | 


on 
= Obes peel — 2) be Cs 


os ih} is 


Y= col 7 =n = const. 


>> > 
ie ae 


Aprés réduction aux variables sans dimension § = 7 v/D —¢ v2/D (abscisse relative 
a l'interphase) et t = / v2/D on peut par une transformation de LAPLACE simplifier 
le probléme et le résoudre. La solution peut étre inversée par une intégration et 
donne une formule assez compliquée pour n(&, t). La concentration de la glace a 
une distance y du point de départ de la cristallisation est alors (avec 4 = y v /D): 


we Sb it AY = 


N,(m) a sn(é= 0,t= 0) 


; = Noo {2 = zerfe|/” — (1— 2s) F = se Erfe (1 — 2s) /2| e-s0-aa) 


aR ‘ F(n, 0) m(t) a, 


esi 
l 4 Fo 
Son Erfex=|/* fer*ay, et 
j —s 


q ple alee ik —s(1—s)n+s¢ n C 
| F(n, 6) = + s) @ > Erfc\ (2s — 1) - + —=——}]. 
| (m, ¢) yx “ ( 2 4 V4 n 


1.2. Le plus fréquemment on part d’une concentration homogéne dans le liquide: 


— an ee 


pee 


| n(x) =O, No FO. 


> eT hee 
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; , 5 oe 
La concentration croit alors de 2,(0) = ms d’abord avec la racine de l’abscisse 


n(n) ~ reo s [1 +.2(1~5) |/ Bap (y < 4) 


puis tend exponentiellement a sa valeur limite 7: 
N,(4) ~ Neo ie (1 — 2s) Qt Rae acct (n > 4). 


Comme la fonction Evfc décroit trés rapidement, et comme en général s est de 
l’ordre du pour cent, on arrive assez vite a la loi exponentielle, et on peut méme 
négliger s par rapport a l’unité. Une fois l’état stationnaire atteint, la glace a la 
méme concentration que le liquide, car il s’est formé a l’interphase une couche de 
concentration S mo, décroissant exponentiellement vers n,, dans le liquide avec la 
largeur caractéristique € = 1, c’est-a-dire ¥ = D/v. 


1.3. Un cas particulier intéressant au point de vue expérimental est celui de la 
solution agitée plus ou moins violement. La diffusion pure se limite alors a une couche — 
d’épaisseur d au voisinage de l’interphase (déterminée par le brassage), la concen- 
tration du liquide étant considérée constante au dela. On peut résoudre facilement — 
le cas stationnaire avec la condition ” = n., pour € = 6 = du/D, obtenant pour la 
concentration dans le liquide a partir de l’interphase : 


n € St; ? = E = ) > 
i tS She coed 
La concentration de la glace est alors 
n(0) Noo Noo 
n= = ; ea 
Ss (S = seee at, Ss 


S* est le coefficient de ségrégation apparent, tendant vers 1 en régime stationnaire 
pour d infini (pas d’agitation) et vers S pour d assez petit (agitation trés violente) 


2. Résultats avec l’acide fluorhydrique 


Les cristaux sont produits dans l'appareil déja décrit [2]. La concentration de 
lacide est déterminée en découpant le cristal cylindrique en tranches de quelques 
millimétres et en mesurant la résistance de l’eau de fonte dans une cellule étalonnée. 
L’influence d’impuretés ioniques accidentelles ne peut pas étre éliminée, et ce n’est 
que pour des concentrations de la glace Supérieures 4 10-7 (mol HF/mol H,O) que 
la dispersion des mesures reste dans des limites acceptables. 


2.1. Sans agitation mécanique on observe une grande influence de la vitesse de 
croissance sur la distribution de l’acide (Figure 1) 


la concentration stationnaire étant d’aprés la thé 
qu’on peut faire passer par les histogrammes concordent au moins qualitativement 


ment des courbes, soit en extrapolant a ¥ 
de 100. Toutefois le dernier de ces procédé 
de croissance trés faible (v < 0,4 b/s) 
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< 
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_ plus grande, car la convection dans le liquide a plus d’influence et suit toutes les 
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variations fortuites de la température a la périphérie du moule. De plus pour des 
concentrations voisines et supérieures 4 10-4, la glace se trouble. On voit que les 
conditions expérimentales ne sont pas favorables, et c’est pour éliminer les influ- 
ences incontrolables du mouvement du liquide que nous avons essayé de l’agiter 


-mécaniquement. 


3:2 4/5 
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Figure 1 Figure 2 
Distribution de l’acide fluorhydrique dans la Coefficient de ségrégation effectif de l’acide 
glace avec la vitesse de croissance comme para- fluorhydrique avec ‘agitation violente de la 
métre, sans agitation de lasolution; 7.,=10~°. solution; tog = 2° 107%. 


2.2. Avec un brassage violent, tel qu’il a déja été décrit [1], la glace est claire, 
mais polycristalline a grain fin, et ne convient donc pas aux mesures électriques. 
Cependant, en reportant les valeurs du coefficient de ségrégation effectif S* en 


fonction de la vitesse dans un diagramme semi-logarithmique (Figure 2), on obtient 


' bien la formule 


ae EN SLT 


log(S* — 1) =log(S—1)—6, 6= a 
en admettant que les conditions du brassage (c’est-a-dire d) et la diffusion ne varient 
pas. L’extrapolation a la vitesse nulle livre S = 60. 


2.3. Un brassage faible a donné les résultats les mieux reproductibles. I] est 
effectué par un disque de polyéthyléne tournant au-dessus de l’interphase a une 
vitesse de quelques t/s. Pour empécher un transport de chaleur par convection dans 
toute la solution et limiter ainsi l’épaisseur de la couche agitée 4 quelques mm, 
le disque est surmonté par un bouchon fixe de polyéthyléne de 2cm de long et 
coulissant avec un jeu de 1 mm dans le moule. Par ce procédé, on élimine la forma- 
tion de bulles (qui apparaissent facilement aux vitesses supérieures a 0,5 u/s sans 
agitation mécanique) mais on limite généralement le volume des unicristaux au 
em’. La distance nécessaire a l’établissement de 1’état stationnaire dépassant 
généralement la longueur du cristal (8-12 cm), on a ainsi une faible variation de 
la concentration de la glace, pratiquement indépendante de la vitesse de croissance 
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(Figure 3) permettant ainsi la détermination de S (Figure 4). On observe une 


variation sensible du coefficient avec la concentration, S montant de la valeur 


50 (Moo < 1075) ala valeur 500 pour x, = 3 - 10-* et descendant ensuite rapidement 
a Dod sour No = 10-*. Dans un diagramme d’état H,O — HF, cette chute correspon- 
drait 4 un rapprochement des deux courbes du solide et du liquide (qui peuvent se 


10-6 


Log S 
4 ; 
a 
HF 
2 
tone bette 
NHg F 
I 
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0,4 u/s 
Log n 
O 
(0) 2) 10 cm => —4 -3 -2 -l 
Figure 3 Figure 4 
Distribution de l’acide fluorhydrique dans la Coefficients de ségrégation en fonction de la 
glace pour deux vitesses de croissance, avec concentration de la solution. (@ = valeur ob- 
agitation faible de la solution; feo = 1074. tenue par DEcROLy et JAccarp [1)). 


joindre 4 un point eutectique). Ce comportement semble donc caractéristique de la 
solution, et laglomeration des impuretés en essaims ou sur des surfaces internes ne 
doit pas y jouer de réle sensible, le trouble observé avec les grandes concentrations 
pouvant trés bien provenir de gouttelettes liquides concentrées prises dans le réseau. 


3. Résultats avec Vammoniaque et le fluorure d’ammonium 


Le coefficient de ségrégation de l'ammoniaque, déterminé par le méme procédé, 
se comporte de la méme fagon (Figure 4), mais atteint des valeurs bien supérieures, 
Pour le fluorure d’ammonium, par contre, il est relativement petit (S < 50) et 
décroit de fagon monotone. Ces résultats peuvent étre comparés a ceux de BRILL 
et ZAROMB [4] qui, pour de faibles concentrations de NH, trouvent de la glace 
pratiquement pure (puisque S ~ 1000) et pour une concentration initiale de 
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2,5: 10-8 trouvent dans de la glace polycristalline 6 - 10-4, ce qui donne S = 42. 
oe ils estiment que les impuretés sont accumulées sur les surfaces, internes, 
enais comme la glace obtenue dans notre cas est presque monocristalline, cette hypo- 
‘thése doit étre rejetée. Il en est de méme de la formation d’essaims pris dans le 
ea car il est peu vraisemblable qu’elle puisse produire une variation si brusque 
mais continue de la ségrégation. Ainsi l’ammoniaque peut donc étre pris en solution 
Bence dans la glace jusqu’a une certaine limite. De plus, aux faibles concentrations, 
_ ilsemble favoriser la formation de glace pure, car des mesures électriques de cristaux 
_ faiblement dotés de NH, et produits rapidement donnent des résultats voisins de 
~ ceux obtenus avec de l’eau distillée, purifiée par échangeur d’ions et congelée trés 
lentement; d’autre part, des bulles de gaz prises dans le cristal ne se produisent 
-méme pas a 3 p/s. 
Les valeurs de S obtenues pour le fluorure d’ammonium sont en général plus 
_ petites que celles de BRILL et ZAROMB, mais si on néglige la correction faite par 
_ ceux-ci et attribuant une certaine proportion de l’azote pris dans la glace a de 
Yammoniaque accumulé sur des surfaces internes, les résultats concordent assez 
_ bien et confirment la grande solubilité du NH,F. 


a2 J 


= 


4. Conclusion 


Il est donc possible de produire des cristaux de glace de quelques cm de longueur 

avec une concentration voulue m,, en utilisant les courbes de la Figure 4 et en 
_ partant d’une solution de concentration 7. telle que Noo|S(Moo) = N. Ll est utile 
- de controler l’agitation du liquide et la vitesse de croissance aussi exactement que 
- possible, en maintenant celle-ci a une valeur inférieure a 0,5 u/s. La facilité avec 
- laquelle les impuretés peuvent étre prises dans le réseau décroit dans l’ordre 
_ NH,F-HF-NH,, et la concentration limite de la solution (pour obtenir des cristaux 
clairs) est dans les trois cas 3 - 10~* environ, sans étre déterminée exactement car 
lappartition des bulles est progressive. 
a Il est remarquable qu’une faible concentration d’ammoniaque permette d’obtenir 
- facilement des cristaux trés purs avec une faible conductivité électrique. Comme 
- acide carbonique est une des impuretés principales de l’eau en contact avec l’air 
- atmosphérique, l’ammoniaque doit certainement jouer un role important dans 
inclusion du COz~. WorKMAN et Reynotps [5], mesurant le potentiel électrique 
durant la croissance, trouvent en effet pour des solutions de NH,, (NH,),COs, 
_ H,CO,, HF et NH,F des tensions de l’eau par rapport a la glace de — 230, — 109, 
0, 0 et + 9 volts respectivement. Il semble ainsi que le ion NH}, bien qu’il puisse 
- a priori entrer facilement dans le réseau, empéche linclusion des anions (OH et 
- CO)z~— et que la différence de potentiel qui en résulte cause a l’interphase un ap- 
_ pauvrissement en cations, et par conséquent un grand coefficient de ségrégation. 
Ceci n’explique pas toutefois le phénoméne pour le HF et le NH,F; la couche limite 
_ pourrait étre partiellement court-cireuitée, diminuant la tension et cet effet pourrait 
- étre également en relation avec l’augmentation et la chute rapide du coefficient de 
_ ségrégation. Seule une étude plus approfondie des propriétés structurelles et ther- 
modynamiques de la surface peuvent donner de plus amples renseignements a ce 
sujet. 
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Summary 


The theoretical distribution of impurities along a growing ice crystal is calculated 
for an arbitrary initial distribution in the liquid, and agrees qualitatively with 
results obtained with hydrofluoric acid. Stirring the liquid phase allows a satis- 
factory determination of the segregation coefficient for HF, NH;,, and NH,F as a 
function of the concentration of the solution. 


(Regu: le 6 septembre 1960.) 


On the Temperature and Humidity Dependence of the Ice Forming — 
Activity of Silver Iodide!) : 


By Urricu Katz, Ziirich?) 


In a recent publication [4]8), it was shown that the ice forming nucleability of a 
substance should be determined as a function of the temperature instead of mea- 
suring only the threshold temperature of the activity. At the same time the author 
described the apparatus and the method by which the activity*) of the copper 
oxides and sulphides had been measured. The interesting results which were then 
reported naturally raised the question of the mode of action of AgI, the best known 
substance for forming ice embryos. Indeed, certain data were already available on 
this subject from other authors: for instance SMITH and HEFFERNAN [9] and 
SOULAGE [10] gave as a function of the temperature, the number of ice crystals 
which are formed in a supercooled cloud per gramm of AgI smoke in the range 
— 4 to — 20°C. The particle size spectra (where these have been measured at all, cf. 
Sano and Fuxuta [8]) are so broad, however, that the number of particles formed 
from one gramm of AgI cannot be determined with sufficient accuracy to enable 
the activity (in our sense) to be calculated. Nevertheless, the opinion was put 
forward that as the temperature is lowered the activity increases exponentially, 
at least above — 10°C, by a factor of about 10 per 2 degrees. 

Our present investigations were carried out with commercial AgI powder with 
which the particles entering the supercooled cloud had the range of sizes shown in 


1) Paper presented at the International Congress of Cloud Physics in Verona, August 9-18, 
1960. 

*) Laboratory of Atmospheric Physics, Swiss Federal Institute of T echnology. 

’) Numbers in brackets refer to References, page 79. 

*) We define the activity as the ratio of the number of ice crystals observed in a supercooled 
cloud to the total number of nucleating particles introduced into the cloud. 
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Figure 1. The measurements, carried out in the same way as described in a previous 
_ paper [4], gave the activity curve shown in Figure 2. The results mentioned above 
_for CuS and Cu,O [4] have been shown in broken lines for the sake of comparison. 
It may be pointed out that in the range —8 to —12°C CuS is just as active as AgI. 

In agreement with earlier authors [3, 7,11] a threshold temperature of — 3 to 
-— 4°C was obtained. These curves (Figure 2) are similar in appearance to those 
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Figure 1 Figure 2 

Relative frequency (on an arbitrary scale) of Ice nucleating activities g of AgI, CuS and 

AglI powder particles of diameter D microns. Cu,O as functions of temperature. 


reported by SouLaGE [10], SmitH and HEFFERNAN [9] and others. The latter results 
have been interpreted theoretically by FLETCHER [2] in terms of the particle size 
distribution of the smokes: below a particle size of approximately 0-5 mw, the thres- 
hold temperature decreases ever more sharply with decreasing particle diameter. 
For smokes with a broad particle size distribution ranging from 1/100 to 1 p, this 
must lead to an effect on the activity such as that observed experimentally. How- 
ever, as we used, deliberately, only particles bigger than 0:5 « for our measure- 
ments, as shown in Figure 1, FLETCHER’S interpretation can only be of minor 
importance in our case: otherwise the activity would have to rise very steeply from 
the threshold temperature in order to reach the maximal value at about — 6°C. 

The small values-of activity which were reached even at relatively low temperat- 
ures are surely the result of a small concentration of suitable sites on the surface 
for the formation of ice embryos. This would be in agreement with the microphoto- 


Saturation with respect to ice 
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graphs shown earlier by Montmory [6] and now also by Bryant, HALLeEt, and 
Mason [1], where the surface concentration of such ice embryo centres on flat AgI 
surfaces is also rather low (order of magnitude about 1073/2). 

Whereas our measurements were, hitherto, always carried out in a supercooled — 
water cloud, i. e. at at least 100% relative humidity, the question of the necessary 
water vapour pressure to form ice embryos was often considered. Assuming that the 
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Vapour pressure of aqueous NaCl solutions; Ice nucleating activities of AgI in clouds of 
numbers on curves indicate g of NaCl in 1000 g aqueous NaCl solutions as functions of tem- 
of water. perature; numbers on curves indicate g of 


NaCl in 1000 g of water. 


ice forming nuclei were only effective where there was water saturation, then ice 
nucleation must become impossible in certain maritime mists, where the drops had 
grown on salt nuclei, and in similar circumstances. In this connection, we carried 
out the following experiments: 

Instead of filling the cloud chamber with a pure water cloud, we used NaCl 
solutions of three different concentrations corresponding to 2, 3 and 5 yu cloud 
droplets grown on 1 w NaCl nuclei. F igure 3 shows the saturation vapour pressures, 
represented as percentages of ice saturation, as functions of temperature for the 
three solutions as well as for pure water. Figure 4 gives the results of the activity 
measurements. The curve for pure water is the same as in Figure 2. For clarity 
the experimental points have not been included but the degree of scatter is the 
same as in Figure 2 except where the curves are broken to indicate a somewhat 


: 
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_ greater uncertainty. This uncertainty may be partly due to the fact that the drops 
increase their concentration as a result of evaporation and thereby lower the vapour 
_ pressure. 

It was naturally to be expected that there would be a clear shift of the threshold 
values to temperatures below the limit of ice saturation; yet the experiment plainly 
shows that the activity is also decreased over the whole temperature range even 

_ for relatively great dilutions. This may be interpreted as meaning that the nucleating 
activity is dependent on the vapour pressure in such a way that activity begins at 
a low percentage of ice supersaturation and increases with rising relative humidity. 

It is interesting to make a comparison with the data on AgI of Mason and VAN DEN 
HEUVEL, according to which water saturation is necessary down to — 12°C, but 

- below that at least 12% supersaturation with respect to ice is necessary; whereas 
we have found activity above — 12°C with slight water undersaturation, and below 
— 12°C with about 5 to 8% ice supersaturation. Our data indicate that AglI acts 
in the whole range of its ice forming activity as a sublimation nucleus. 

These results are still too sparse and not accurate enough to form the basis of 
theoretical calculations about the influence of vapour pressure, but our experiment, 
which was only designed to furnish some general indications for a future line of 
research, does show how relative humidity plays a part in the question of icing 
nucleability and that, in certain cases, this matter should be taken into considera- 
tion in meteorological applications. 
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Zusammenfassung 


An AglI-Pulver mit der Korngréssen-Verteilung nach Figur 1 wurde die Eis- 
keimbildungsaktivitat in Funktion der Temperatur gemessen ; Figur 2 zeigt das 
Resultat, das kurz diskutiert und mit Experiment und Theorie anderer Autoren 
verglichen wird. 

Die Abhangigkeit der Aktivitat vom Dampfdruck zeigte sich in einem orien- 
tierenden Versuch, bei dem der unterkiihlte Nebel aus Steinsalzlosungstrépfchen 
bestand; mit den so reduzierten Dampfdriicken gemass Figur 3 ergaben sich die 
in Figur 4 dargestellten Aktivitatskurven, woraus zu schliessen ist, dass AgI im 
ganzen Bereich seiner Eiskeimbildungsaktivitat als Sublimationskern wirkt. 


(Received: August 12, 1960.) 
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Heat Transfer from Laminar Flow Through Cylindrical Tubes 


By SERGEJ Panor and JANEZ STRNAD, Ljubljana‘), Yugoslavia 


The differential equation for heat transfer in laminar flow is 


is 
divgrad T = © 2 (5 + uv grad T) (1) 


VO tin tetas inl 


with all physical properties (g density, c, specific heat, 2 thermal conductivity) of 
the fluid assumed constant. For stationary flow through a cylindrical tube (dia-— 


meter 27) with constant wall temperature T,, the solution is sought in the form 


T= Ly LCs R,(4) exp (— by 2) (2) . 


where w is ¥ = r/r). C, are constants depending on the initial condition at the — 


entrance of the tube. The equation for R,,(x) is 


if 4 
Rela) + Ral) + (fie + 3 (1 —24)) Ryle) = 0 (3) 
with the boundary conditions 
RNa RH =O. (4) 


Pe stands for the Peclet number Pe = 2% 7, Cy g/A (u is the mean velocity) and p2_ 


for the eigenvalue p2 = b, 7, Pe. At a sufficient distance from the entrance of the 
tube the temperature profile can be approximated with the first eigenfunction 
R,(*) = R(x). The first eigenvalue is p? = 2. In most cases of technical interest 
the term p4/Pe? can be neglected. The solutions of equation (3) and of the correspond- 
ing equation for p4/Pe? = 0 have been calculated by various methods to different 
degrees of approximation [1] — [4]?), 

The solutions of (3) can also be expressed with confluent hypergeometric 
functions. This is of practical interest, becausethe confluent hypergeometric 


functions are quite well tabulated [5]. The solution R(*) expressed with the confluent 
hypergeometric function is 


where 


2° 


The first boundary condition is satisfied already and from the second boundary 
condition it follows 


om $ (1+ 4) -4 


(= @; 1) p) == (5) 


Using tables [5] the calculation of the relation between the first eigenvalue p2 and 
Pe is straightforward. 


1) Institute of Physics, University of Ljubljana. 
*) Numbers in brackets refer to references at the end of the paper. 
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The Nusselt number is defined [3, 4] by 


The value of K 


1 at Sil 
(Pe) | Ra) nae R(x) (1 — 2?) x a) 
Naseeal 


is almost constant, e.g. K(0) = 1-445 and K(oo) = 1-419. Therefore we get quite 
a good approximation for Nu, at least for great values of Pe, if we insert in the 
_ equation (6) K = K(oo). If greater accuracy is desired the value of Nu for various 
_ Pecan be calculated with the aid of tables. 


: 


Nu(0)=4-1805 | 4.2 


Figure 
The functions p2(Pe) and Nu(Pe). 
~ Nu(Pe) was calculated from equation (6) with an approximation K(Pe) = K(oo) (1 + 0-0032 p4/Pe?). 
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Zusammenfassung 


Die Differentialgleichung fiir Warmeleitung in einer durch ein zylindrisches 
Rohr strémenden Fliissigkeit, die wie bekannt schon 6fters behandelt wurde, wird 
in der vorliegenden Arbeit mittels der konfluenten hypergeometrischen Funktion 

' gelst. Bei der Berechnung des ersten Eigenwertes wird von den neulich heraus- 
gegebenen Tafeln Gebrauch gemacht. 


. (Received, August 8, 1960.) 
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Electron Density Distribution in a Cylinder’) 


By Hetmut D. Weymann, College Park, Md., USA?) 


Introduction 


In experimental studies of electron diffusion from argon plasmas in shock tubes 
it was found that appreciable electron concentrations can be measured up to 1m 
ahead of the shock front [1]*). A complete theoretical description of this nonstation- 


ary diffusion process is very difficult, as the diffusion differential equation is non- 


linear due to the electric force on the electrons. The electric field, generated by the 
charge separation, is moreover, strongly three-dimensional. There is, however, one 
aspect of the problem which is amenable to a theoretical treatment. As the tube 
diameter is very small (about 3cm) compared with the axial extension of the 
phenomenon (about 100 cm) it may be assumed that the diffusing electrons will 
establish a quasi-equilibrium density distribution in the radial direction in all but 
the region directly ahead of the shock. The radial distribution at a given point and 
instant can then be approximated by the stationary radial distribution of electrons 
in an infinitely long cylinder provided the mean density over the cross section is 
the same in both cases. This is the problem considered here. 


Electron Density Distribution 


The aim of the following calculations is to determine the equilibrium density | 


distribution of electrons of a given temperature T in an infinitely long insulating 
cylindrical tube of radius FR, which is capable of withstanding any field strength. 
There may be a background of neutral gas in the cylinder. It is, however, assumed 
that no positive ions are present. 


The differential equation for our problem can be found easily. It is known that 
the electron density as a function of the (diffusion) potential is given by a Boltzman 
distribution*) 


e 
n, = n* exp (F) : (1) 


(y = potential, # = Boltzmann constant, T = electron temperature, e = charge of 
an electron, n* = electron density at y = 0). In addition one has Porsson’s equation 
which, in the absence of positive ions, assumes the form 


en, 
AN) yo (2) 
2 eo 
1) This research was supported in part by the United States Air Force through the Air Force 


Office of Scientific Research of the Air Research and Development Command under Contract 
No. AF 49(638)-401, 

) Institute for Fluid Dynamics and Applied Mathematics, University of Maryland; now at 
the University of Rochester, Rochester, N. Ye, WSL 


8) Numbers in brackets refer to References, page 85. 
4 eaters aaa 
) In equilibrium the pressure force on the electrons: — grad p = —k T grad n is equal and 


opposite to the electrical force — ne E = ne grad y. Hence n~ grad n = (e/kT) grad yw, which, for 
T = const, can be integrated to give equation (1). 


j 
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) (The dielectric constant of the background gas has been set equal to unity.) The 
_ Laplacian operator in cylindrical coordinates is 


y oe 1 0 0° Nati 
- A = <—— —— ee ee : 
oat FoF dg? T Og? To 9° 
_ With the assumed cylindrical symmetry and independence of z equations (1) and 
'/(2) then result in 
O?y 1 Op en* ew 
Ay = : Set 
| EPR ag chigo™ 55 exp (5 r) i ie 
_ By substituting 
Oe ok e2 n* bos 
i ie and Ser A (4) 
into equation (3) we obtain 
i eat ae (5) 


Equations of this type have been investigated by LEMKr [2]. He showed that equa- 
_ tion (5) can be brought into the form 


ib uw G 
Ua = IN. ees ete 10) 6 
uU 3 (w’) 5 (6) 


(C = constant of integration not directly related to A) 
which is a Riccati differential equation in wu’. With the transformation 


u=—2\inw+Int (t= constant) (7) 


equation (6) goes over into the linear differential equation 


w’ we 
" re 8 
w 5 + 2 0. (8) 
A solution of equation (8) is 
w= 7,0 + 720" (9) 


where f, » are the solutions of the characteristic equation 


: oy 
pga Os a = 0 


and hence 
f 


By 2 = 1 = (2 = F : (10) 


The value of 7 in equation (7) can be determined by substituting equation (9) 
into equation (7) and in turn equation (7) into equation (5). This results in 


4—2C 
i. aig ad Neonat, Terai (11) 
Hence from equations (7), (9), and (11) 
j 4—2€° 
“= — 21n (v1 + 72.0%) + In(— Daa ). (12) 


The value of the integration constant C can now be determined by the following 
consideration. As the potential y = uk T/e must remain finite at g = 0 either f, 


= Y 


’ 
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or B, has to be zero and the other must be greater than zero (see equation (12)). A 
comparison with equation (10) shows that the only value of C for which one of the 
is zero is C = 0. Then f, = 0 and f, = 2. For convenience we shall assume further 
that y = u = 0 for @ = 0. From equation (12) follows then - 


y= — ae (13a) 
and with equation (4) 
ial } eat 
A= aa (13b) 
Vo + 8e,kT 


Introducing these relations for f,, B,, C, and y,/y, into equation (12) gives 


eh Le ae eae 
2In (1 at (14a) 
where 
re ae he 
B= ae (14b) 


the Debye length, is defined with n* the electron density for y = 0 at 9 = 0. 


The electron density distribution is found by substituting equation (14a) into 


equation (1) which gives 
n 


We (1 — 92/8 is 2, 


(15) 


where n* and /p are quantities which as yet have to be determined. As , has to 
remain finite for 0 << g < R, the Debye length can assume only values which 
satisfy the inequality 


82 <eilce (16) 
or with equation (14b) 
Ae 8 €o RD * 
ie Wow aa max (17) 


athe maximum value which the electron density can assume at 9 = 0 is, therefore, 
determined by the temperature and the tube radius. 


Now we shall derive a relation between n* and the mean density for a cross 
section 
R 
. y 1 
n 7 R ‘ 220 do n,(o) . (18) 


0 


After introducing equation (15) into equation (18) and integrating, one obtains 


n* 
n= = R2/8 i, a (19) 
By using equation (14b). one finally arrives at 


1+ Re2n/BkT e, ’ 


(20) 


where n* the electron density at @ = 0 is given as a function of the mean density n, 
the cylinder radius R and the electron temperature T. 


3 
i 
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In Figure 1 several density profiles have been plotted for a cylinder of radius 
1cm and an electron temperature of T = 104° K. For-these values /%, has to be 
_ larger than 0-125 cm*, which in turn means that the upper limit for n* is 
nx yy = 3-81 X 108 cm. 

_ Figure 1 shows that for m < n*_ the electron density is nearly independent of g. 
For n > n*_, the density at the center approaches n* , while the density at the 
wall increases with the square of m which can be verified by substituting equations 

_ (14b) and (20) into equation (15). This results in 


R2 e2 *) n2 PR? e2 
wy 


8h Te, Snr eae (21) 


L=1001cm 
n-19x10° cm=? 


L=11¢m 
N=1.81%10’ cm 
=1.20M 
= 8.63>10 cm * 


L=2cm 
n=127«10%cm-* 


L=5em 
N=158~10°cm=$ 


Figure | 


Radial electron density profiles in an infinitely long cylinder (Electron temperature T,= 104° K; 
L?=8 7%). 
D 
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Zusammenfassung 


Die radiale Elektronendichteverteilung in einem unendlich langen Zylinder wird 
fiir den Gleichgewichtsfall berechnet. Dabei ergibt sich die Existenz einer maxi- 


malen Dichte n* welchem Wert die Dichte in der Achse des Zylinders zustrebt, 


wenn die mittlere Dichte 7 im Querschnitt gegen unendlich geht. Figur 1 zeigt, 


dass fiir » < m*_ die Dichte iiber den Querschnitt praktisch konstant ist. Fir 


n > nx, Steigt die Dichte an der Wand mit 7? an. 


(Received: July 5, 1960). 
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Symposium on Electromagnetics and Fluid Dynamics 
of Gaseous Plasma in New York City 


The eleventh annual international Polytechnic Institute of Brooklyn Sympo- 
sium, scheduled for April 4-6, 1961 in New York City, will be devoted to ‘Electro- 
magnetics and Fluid Dynamics of Gaseous Plasma’. The area of interest covers 
collective electromagnetic and fluid dynamic phenomena and their self-consistent 
analytical description in the temperature range characteristic of the transition 
between partially and fully ionized gaseous plasmas. 

The symposium is being organized by the Department of Aerospace Engineering 
and the Microwave Research Institute of the Polytechnic Institute of Brooklyn in 
cooperation with the Institute of Radio Engineers and the Institute of the Aero- 
nautical Sciences and with the co-sponsorship of the Air Force Office of Scientific 
Research, the Office of Naval Research, and the U. S. Army Signal Corps. 

Symposium Committee 
Polytechnic Institute of Brooklyn 
55 Johnson Street, 

Brooklyn 1, New York 


Mathematisches Symposium des internationalen Rechenzentrums 
der UNESCO (PICC) in Rom, 22. bis 24. September 1960 


__ Das Thema des Symposiums war die numerische Behandlung von gewéhn- 
lichen Differentialgleichungen und von Integral- und Integro-Differentialglei- 
chungen. Hauptvortrage wurden von Prof. R. CourANtT (New York), Herrn 
Genuys (Paris) und Prof. WaLTHER (Darmstadt) gehalten und ausserdem etwa 
50 spezielle Vortrage von Fachleuten aus den verschiedensten Landern. 


Der Kongressbericht (Proceedings) wird demndchst im Birkhauser Verlag 
(Basel und Stuttgart) erscheinen, 


: 


s 
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Buchbesprechungen — Book Reviews — Notices bibliographiques 


Mathematics in Physics and Engineering. Von J. Irvine und N. Mut- 
LINEUX (Academic Press, New York und London 1959). 883 S., 105 Abb.; $11.50. 
_ Als Band 6 der Reihe «Pure and Applied Physics» erschienen, méchte das Buch 
von IrRvinG und MULLINEUX dem angehenden Physiker und Ingenieur das heut- 
zutage erforderliche Riistzeug an Kenntnissen aus der Mathematik vermitteln. Das 

- Hauptgewicht der Darstellung wird von den Verfassern naturgemadss auf die An- 
wendung der verschiedenen mathematischen Methoden auf die Probleme der 
_ Physik gelegt, wobei diese dem Gebiet der Elastizitat, der Uberschallstrémung, des 
Elektromagnetismus, der Wellenmechanik und der Warmeleitung entnommen 
wurden, und weniger auf rein mathematische Fragen. Einige weitergehende mathe- 
matische Ausfiihrungen wurden deshalb in einen Anhang verwiesen. Im wesent- 
lichen sind die folgenden Gebiete behandelt: Gewohnliche und partielle Differential- 
gleichungen, spezielle Funktionen, Transformationstheorie, Matrizen, Variations- 
technung, Funktionentheorie, Differenzenrechnung, Integralgleichungen und nu- 
merische Methoden. Diese Ubersicht zeigt, dass der Rahmen recht weit gespannt 
ist, aber dem entspricht, was heute gefordert werden muss. Als Lehrbuch enthalt 
das Werk zu allen Abschnitten neben vielen, mehr oder weniger weitgehend durch- 
gerechneten Beispielen auch zahlreiche Ubungsaufgaben (deren Losungen anhangs- 
weise beigefiigt sind) sowie einige ausgewahlte bibliographische Hinweise auf aus- 
fiihrlichere Spezialwerke. — Als Ganzes gesehen, handelt es sich um ein recht sorg- 
faltig und auf Grund einer reichen Erfahrung zusammengestelltes einfiihrendes 


Werk, das den angestrebten Zweck sehr wohl zu erfiillen vermag. 
E, RotH-DESMEULES 


Thermische Turbomaschinen. Dampfturbinen, Gasturbinen, Turboverdich- 
ter. — Von WALTER TRAUPEL. Bd. 1: Thermodynamisch-strOmungstechnische Be- 
rechnung (Springer-Verlag, Berlin 1958). 407 Seiten, 402 Abb., sFr. 62.85. 

Prof. Dr. WALTER TRAUPEL, Inhaber des Lehrstuhls fiir thermischen Maschi- 
nenbau an der ETH in Ziirich, hat unter Verarbeitung letzter Erkenntnisse aus 
Thermodynamik und Strémungslehre ein neues Standardwerk tiber den thermi- 
schen Turbomaschinenbau herausgegeben. Im ersten Band werden in elf Kapiteln, 
welche in sich abgeschlossene Einheiten bilden, die thermodynamischen Grundlagen, 
die Berechnung der Arbeitsprozesse, die stromungstechnischen Grundlagen, die 
Arbeitsverfahren thermischer Turbomaschinen, die elementare Theorie der Stufe, 
das Schaufelgitter, die réumliche Stromung durch Turbomaschinen, die Berech- 
nungsunterlagen, die Auslegung der mehrstufigen Maschine, Wellendichtungen und 
Schubausgleich sowie das Verhalten der Turbomaschinen unter geanderten Be- 
triebsbedingungen behandelt. Jedem Abschnitt ist ein Literaturverzeichnis bei- 
gefiigt, welches dem Fachmann ein vertieftes Studium der in ihrer Gesamtheit darge- 
stellten Probleme erméglicht. Von unmittelbarer Aktualitat sind, im Zusammen- 
hang mit dem Einsatz der Atomenergie, die Ausfiihrungen iiber Nassdampfturbinen. 
Turbine und Verdichter werden auf einheitlicher Grundlage véllig analog behandelt, 
was gegeniiber den bekannten historischen Darstellungsarten einige Umstellungen 
_ bringt. Das Buch von Traupel zeichnet sich durch einen streng logischen Aufbau, 

seine klare Sprache und die systematische Anwendung des MKSA-Systems aus. 
Finleitend werden die einfachen Regeln fiir die Umrechnung in das technische 
“MaBsystem mitgeteilt. Die erwahnten Vorziige erleichtern dem Studierenden, wel- 
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cher sich erstenal at dem thermischen Turbomaschinenbau zu befassen has, das 
Eindringen in’ die neue Materie ausserordentlich. Aber auch dem in der Industrie | 
tatigen Ingenieur bietet das Buch von Traupel durch seine exakte Erfassung der 
Einzelprobleme eine Fiille von Anregungen. Die praktische Anwendung der ent- 
wickelten verfeinerten Theorie diirfte zu weiteren Verbesserungen der thermischen 
Maschinen fiihren. In einer Zeit, die ein Uberangebot an beschreibenden Sammel- 
werken kennt, bedeutet das Lehrbuch von Traupel, welches in gedrangter und gut 
verstandlicher Form jeweils das Wesentliche der einzelnen Probleme erfasst, einen : 
Lichtblick. W. R. Duss” 


Praxis der konformen Abbildung. Von W. von Koprenrets und F. STALL- | 
MANN ((Band 100 der Grundlehren der Mathematischen Wissenschaften] Springer- 
Verlag, Berlin, Géttingen und Heidelberg 1959). 375 S., 250 Abb.; DM 69.-. 

Das urspriinglich von W. von KoppENFELS begonnene und in erster Linie als” 
Nachschlagewerk geplante Buch, dessen Bearbeitung schliesslich StaLLMANN be- 
sorgte, vermittelt in einem ersten, sorgfaltig ausgewogenen Teil das Wichtigste aus _ 
der allgemeinen Theorie der konformen Abbildung, namlich iiber komplexe Funk- 
tionen, konforme Abbildungen durch spezielle Funktionen, Potential- und Funk- | 
tionentheorie, Theorie und Praxis der Polygonabbildungén und die verschiedenen 
Naherungsverfahren. Der zweiteTeil enthalt dann den Katalog der wichtigsten kon- 
formen Abbildungen. Dabei werden die Abbildungsfunktionen explizit angegeben 
und, soweit dies iiberhaupt méglich ist, auch die Abhangigkeit der Lésungen von 
den Parametern untersucht. Die Zusammenstellung erstreckt sich auf Zweiecke, — 
Geraden- und Kreisbogendreiecke, ein- und zweiteilige Geradenvierecke, Kreis-_ 
bogenvierecke, Kegelschnittpolygone und -schlitze, Parallel- und Geradenschlitze, _ 
konzentrische Kreisbogenschlitze, Die sorgfaltige und eingehende Darstellung © 
macht das Buch, dessen Ausstattung hervorragend ist, zu einem wichtigen Hand- 
buch und Nachschlagewerk, das fiir die Praxis der konformen Abbildung ein un-— 
entbehrliches Hilfsmittel sein wird. E. RotH-DESMEULES 


Strukturtheorie der Wahrscheinlichkeitsfelder und -raume. Heft 24: | 
Ergebnisse der Mathematik und ihrer Grenzgebiete. Von D. A. Kappos (Springer- 


Der Verfasser hat schon 1948 [Zur Math. Begriindung der Wahrscheinlich- 
keitsrechnung, Sitzungsberichte der Bayrischen Akad. der Wissenschaften, 1948] 
ein Wahrscheinlichkeitsfeld erklart als eine additive, normierte, strikt positive | 
reelle Funktion auf einem Boolring mit Eins, und dabei die Vorteile dieser Metho- _ 
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definiert ist. Die Kapitel 4 bis 7 behandeln kartesische Produkte und verschiedene 


Kapitel noch der Begriff der bedingten 
H.MatziIncEeRr 


apitel tiber topologische bzw. kompakte _ 


Neuerscheinung eee. 
KOWALSKY 


TOPOLOGISCHE RAUME 


Von Prof. Dr. H.-J. Kowalsky, Professor an der Universitat Erlangen 
271 Seiten, Preis Fr./DM 40.— 


Mathematische Reihe, Band 26, Sammlung «Lehrbiicher und Monographien 
aus dem Gebtete der exakten Wissenschaften» 


Inhalt: I. Grundlagen: Mengen, Verbande, Filter. II. Topologische Raume: 
Zusammenhang zwischen Topologie und Grenzwert, Grundbegriffe, Tren- 
nungseigenschaften, Machtigkeitsbedingungen. III. Kompaktheit und Zu- 
sammenhang (einschliesslich lokaler Begriffe und parakompakter Raume). 
IV. Abbildungen: Stetige, offene und abgeschlossene Abbildungen, Homéo- 
morphie, vollstandig regulare Raume, Quotienten-, Produkt-, Summen- und 
Abbildungsrdume. V. Erweiterung und Kennzeichnung topologischer 
Raume: Allgemeines Erweiterungsprinzip, Kompaktifizierung, Einbettungs- 
und Darstellungssdtze. VI. Metrische und uniforme Raume: Metrisierung, 
Gleichmassigkeit, Vervollstandigung. VII. Topologische Gruppen, Anwen- 
dungen: Grundbegriffe aus der Theorie der topologischen Gruppen, Ringe, 
Korper und Vektorraume, Approximationssatz von Stone-Weierstrass, 


induktiver und inverser Limes. 


Das vorliegende Lehrbuch soll in die mengentheoretische Topologie einfiih- 
ren. Ausgangspunkt ist eine Verallgemeinerung des Grenzwertbegriffs, die 
einen natiirlichen Zugang zu den topologischen Raumen vermittelt. Kenn- 
zeichnend fiir die Darstellung ist die weitgehende Verwendung der Filter, 
die nicht nur eine tibersichtliche Schreibweise und Begriffsbildung, sondern 
auch eine formale Beweisfiihrung erméglichen. Die Stoffauswahl wird durch 
Hinweise und Aufgaben am Ende-eines jeden Paragraphen erganzt. 


Zu beziehen durch Ihve Buchhandlung — Obtainable from your bookseller — 
Commandes & votre libraire 


Birkhauser Verlag - Basel und Stuttgart 


Neuerscheinung 


Statistische Methoden 


fiir Naturwissenschafter, Mediziner und Ingenieure 
von 


ARTHUR LINDER 
Professor an der Universitat Genf und an der Eidgendssischen Technischen Hochschule Ziirich 


Dritte, umgearbeitete und stark erweiterte Ausgabe 
a ts ‘ (1960) 484 Seiten mit 58 Figuren und 73 Beispielen. In Ganzleinen Fr. 54.— (DM 54.-) 


Mathematische Reihe, Band 3, Sammlung «Lehrbiicher und Monographien aus dem Gebiete der exakten 
Wissenschaften» 


' _ Angesichts der unvermindert anhaltenden Nachfrage hat sich der Verfasser entschlossen, seine bekannte 

* Einfiihrung in die neueren Methoden der mathematischen Statistik stark zu erweitern. Gleichzeitig wurde 
; . der den Anwendungen gewidmete Teil des Buches véllig umgestaltet. 

Neu aufgenommen wurden unter anderem: Verschiedene Anwendungen yon Chiquadrat, die Streuungs- 

zerlegung bei ungleichen Klassenhaufigkeiten, die Bestimmung von Streuungskomponenten, die nichtlineare 

a - Regression, die «analysis of covariances, das Schatzen von Parametern und einige Transformationen von 

ij |  Prozentzahlen (Arc sin, Probit, Logit, Loglog). 


4 : 1959 erschien vom gleichen Autor: 


Planen und Auswerten von Versuchen 


Eine Einfiihrung fiir Naturwissenschafter, Mediziner und Ingenieure 
Ta? 2. Auflage, 182 Seiten mit 9 Fig., Ganzleinen Fr./DM 21.— (1959) 


| Reihe der experimentellen Biologie Band 13, Sammlung «Lehrbiicher und Monographien aus dem Gebiete 
fe der exakten Wissenschaften» 


f Diese Einfithrung richtet sich an Naturwissenschafter, Mediziner und Ingenieure; sie setzt keine Kenntnis 

: in mathematischer Statistik voraus. Der Leser wird angeleitet, Versuche richtig zu planen und einwandfrei 

auszuwerten, Die grimdlich durchgearbeiteten Anwendungsbeispiele aus der biologischen, medizinischen, 
industriellen und landwirtschaftlichen Forschung bilden einen wichtigen Bestandteil des Buches. 


1961 erscheint vom gleichen Autor: 


Handliche Sammlung 
5 mathematisch-statistischer Tafeln ~ 


40 Seiten, Fr./DM 4.50 


Die vorliegende Sammlung statistischer Tafeln soll die umfangreicheren, bisher erschienen nicht ersetzen: 
sie stellt vielmehr einen handlichen Auszug der am meisten verwendeten Tafeln dar, die dazu bestimmt sind, 
stets griffbereit auf dem Arbeitstisch zu liegen. Ausser den Tafeln von t, x2 und F enthiilt diese Sammlung 
Tafeln fiir die Transformationen von Prozentzahlen mittels Arcus sinus, Probits, Logits und Loglog, sowie 


Tafeln fiir zufallig angeordnete Zahlen, die besonders beim Planen von Versuchen und Stichprobenerhebungen — 
zu beniitzen sind. : 


Zu beziehen durch Ihre Buchhandlung -— Obtainable from your bookseller 


BIRKHAUSER VERLAG - BASEL UND STUTTGART 
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